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PREFACE 


Ballistic  Missile  Defense  impact  lethality  problems  involving  non¬ 
nuclear  kill  mechanisms  and  fragment  damage  processes  have  become 
Increasingly  complex.  Complicated  target  configurations,  extreme  impact 
conditions,  and  varying  angles  of  attack  require  the  use  of  sophisticated 
computer  simulations  to  model  problem  geometries.  Several  different  hydro¬ 
codes  have  been  used,  but  where  complicated  geometry  existed,  modifications 
have  had  to  be  made  to  the  configuration  to  fit  the  constraints  of  the 
codes.  In  order  to  analyse  these  problems  without  the  need  of  bending  the 
problem  to  fit  the  code,  it  was  decided  to  construct  a  new  code  that  would 
handle  these  configurations  without  alteration.  Rather  than  develop  a  new 
code  from  scratch,  an  existing  two-dimensional  code  was  chosen  as  a  founda¬ 
tion,  where  previous  work  and  experience  could  be  built  upon.  Since  in 
previous  work  with  two-dimensional  hydrocodes,  it  was  determined  that  one 
code,  CSQII,  performed  better  than  others  In  correlating  with  test  results, 
it  was  decided  to  use  this  code  as  a  base.  The  result  was  K3,  a  computer 
code  which  simulates  material  response  in  three-dimensions.  This  report 
provides  the  technical  information  required  to  understand  the  CSQII  to 
K3  extension,  details  the  Inputs  and  co  trol  cards  necessary  for  execution 
on  a  CDC-7600,  and  presents  the  results  from  program  development  and 
checkout. 

The  authors  wish  to  extend  special  thanks  to  Mr.  S.L.  Thompson  of 
Sandia  Laboratories  for  his  ideas  and  comments  during  the  development  of 
K3.  His  help  greatly  aided  the  designers  and  shortened  the  K3  development 
period.  Also,  the  authors  are  indebted  to  the  many  people  who  contributed 
to  the  K3  effort  by  providing  help  In  debugging  the  code. 
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I  INTRODUCTION 

K3  is  a  computer  code  developed  by  Kaman  Sciences  Corporation  under 
contract  to  the  Ballistic  Missile  Defense  Advanced  Technology  Center  for  the 
simulation  of  three-dimensional  material  motion.  K3  is  an  extension  of  CSQ, 
a  program  for  computation  of  two-dimensional  material  motion  developed  by 
Sandia  Laboratories,  Albuquerque,  New  Mexico.  The  form  of  CSQ  has  been 
retained  in  the  computational  method,  input  and  output  form  and  in  the  names 
of  subroutines.  A  familiarity  with  CSQ  would  greatly  aid  the  users  of  K3 
since  the  basic  material  explaining  many  concepts  common  to  the  two  codes  is 
not  repeated  here.  Rather,  it  is  the  intent  of  this  report  to  present  where 
and  why  modifications  have  been  made  in  order  to  make  the  extension  to  the 
third  dimension.  Users  should  note  that  the  K3  program  is  not  perfectly  clean 
in  the  sense  that  all  program  references  to  "CSQ"  have  been  converted  to  "K3". 
The  development  of  K3  is  an  ongoing  task  where  improvements  are  continually 
being  made.  In  the  course  of  this  work  it  is  hoped  that  eventually  K 3  will 
cease  to  be  a  hybrid  and  will  become  a  unique  program. 

K3  is  a  title  that  applies  both  to  a  program  and  to  the  package  of  programs 
which  include  and  support  K3.  The  support  programs  are 

PREK3  -  generates  the  COMDECKs  used  by  other  programs  in 

the  package 

K36EN  -  generates  the  K3  (re) start  tape  from  the  initial 

problem  configuration 

K3PLT  -  graphics  program  for  plotting  K3  results 

Other  programs  which  are  part  of  the  CSQ  package  but  not  listed  here 
have  not  yet  been  converted  to  the  K3  system.  Also,  several  capabilities  of 
CSQ  have  yet  to  be  included  within  the  K3  programs.  These  include 

.  cylindrical  geometry 
.  rezone/ remeshing  operations 
.  energy  sources 
.  high  explosives 


t 


.  tracer  particles 
.  radiation 

.  various  plot  options  in  K3PLT 

Wherever  possible,  these  CSQ  capabilities  which  are  embodied  in  specific 
subroutines  have  been  retained  pending  future  conversion  and  modification  efforts. 
Since  these  routines  still  exist  in  the  code,  it  is  theoretically  possible  to 
execute  these  capabilities,  but  since  results  are  neither  known  nor  guaranteed, 
use  of  these  capabilities  is  highly  discouraged.  Wherever  certain  inputs 
could  route  the  execution  path  through  unmodified  portions  of  the  code,  comments 
are  included  to  warn  the  user  that  only  the  K3  standard  input  values  should  be 
defined. 

This  report  is  divided  into  these  main  parts:  Part  I,  Introduction;  Part  II, 
Sections  A  and  B,  is  intended  to  aid  those  users  interested  in  the  physical 
considerations  inherent  in  the  K3  code.  Section  A  presents,  without  derivation, 
some  of  the  more  important  modified  and  generalized  equations  which  form  the 
basis  of  the  K3  program.  Section  B  attempts  to  give  the  user  an  understanding 
of  how  the  physical  laws  governing  material  flow  are  incorporated  into  K3 
through  use  of  finite  difference  relationships  and  the  spatial  mesh.  Though  it 
is  not  necessary  for  the  user  to  understand  the  equations  presented  or  the  code 
structure  to  successfully  use  the  K3  program,  a  review  of  the  material  of  these 
sections  should  be  made  for  better  understanding  the  necessary  inputs  and  the 
resultant  outputs.  The  next  part  is  a  programmers  guide  presenting  a  description 
of  the  restart  tape  and  COMMON  blocks  for  use  by  users  who  wish  to  modify  the 
program. 

Section  IV  is  a  user's  guide,  with  a  listing  of  control  cards  necessary 
for  code  execution  (Section  A),  and  inputs  for  the  programs  in  the  K3  package 
(Section  B).  Appendix  A  gives  a  brief  overview  of  the  procedures  used  and 
results  obtained  in  the  course  of  K3  program  development,  test  and  debug. 
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II  PHYSICAL  AND  NUMERICAL  RELATIONSHIPS 
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A.  PHYSICAL  RELATIONSHIPS 
1 .  Conservation  Equations 


The  Lagrangian  forms  of  the  conservation  taws  in  Cartesian  coordinates 
which  control  the  three  dimensional  flow  of  substances  in  K3  are  as  follows: 

Conservation  of  Mass 


1  aP  3VX 
p  It  -  3x 


(2.1) 


Conservation  of  Momentum 


3Vy  a 
-P — -  =  - 


3a 


=  -4r(a  +p  +  Q  )  +  + 

ax  v  xx  r  wxx'  ay  az 


3o 


xz 


(2.2) 


[ 

[ 

l 

I 

t 

! 

K 

I 

I 


3V  3a, 
■p  It  =  lx 


3x  +  3y  ^ayy  +  Pr  +  ^yy^  + 


3a 

lz 


(2.3) 


h  (a 


zz 


Pr  + 


<>zz> 


(2.4) 


Conservation  of  Energy 
E_ 


It  <Em  + 


-<Pm  +  pr  +  5)  h 


i  V.  T  +  r  +  pd 


(2.5) 


In  the  preceding  conservation  equations  t  denotes  time  and  p  denotes  density. 
Assuming  that  the  index  notation  i  *  1,  2,  3  corresponds  to  Cartesian  coordinates 
x,  y,  z,  then  V^,  a^  and  Q^.  denote  the  velocity,  stress  tensor  and  artificial 
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viscosity  tensor  components,  respectively.  The  notation  pertaining  to  the 
remaining  variables  is  as  follows: 

Em  -  specific  material  internal  energy 

Er  -  radiation  energy  density 

Pd  -  stress  deviator  power 

Pm  -  material  pressure 

Pp  -  radiation  pressure 

Q  -  artificial  viscosity  quantity  (Section  II-A-2) 

F  -  energy  flow  flux  vector 

• 

f  -  specific  internal  energy  production  rate 

The  three  dimensional  conservation  laws  can  be  rewritten  like  those 
pertaining  to  CSQ.  The  procedure  for  rewriting  these  equations  is  contained 
in  Section  II-l,  Reference  1. 


2.  Artificial  Viscosity 

The  conservation  of  mass  equation  (2.1)  can  be  written  as 


I  3p  _ 
P  "3t  “ 


V.  V  = 


f  +  f  +  f 
x  y  z 


where 


f 


x 


avu 
?  =  l 
y  3y 


(2.6) 


(2.7) 
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The  generalization  of  the  CSQ  viscosity  tensor  to  three  dimensions  using 


I 

(2.7)  is  as  follows 

.  Define 

I 

fx  ■  ~fx 

if  fx  <  0 

I 

=  0 

if  fx  >  0 

(2.8) 

r 

II 

** 

if  fy  <  0 

i 

=  0 

if~fy>0 

(2.9) 

i 

fz  "  ~fz 

if  fz  <  0 

i 

=  0 

if  fz  >  0 

(2.10) 

With  these  functions  the  diagonal  elements  of  the  artificial  viscosity  tensor 
'  are  nonzero  when  the  material  is  undergoing  compression  and  are  given  by 

«x*  '  ♦  fy  +  fz>  •  B«cs]fx  <2-"> 

T* 

1  qyy‘‘‘[B5(fx+ Vfz>  -  Vs]fy  <2->2> 

r 

I*  Q„  ■  »[B2(fx  *  fy  *  fz)  -  BtCs]f2  (2.13) 

where  B.  and  B„  are  constants. 

[  q 

f 

Q 

B 


-7- 


The  artificial  viscosity  tensor  components  Qxy,  Qyx,  Qxz,  Qzx,  Qyz  and 
Qzy  are  all  zero  in  the  K3  program.  The  quantity  Q  in  equation  (2.5)  is 


0  f  + 

_  WXX  X 


Q  f 


+  Q  f 
yzz  z 


f  +  f 
X  y 


+  f. 


(2.14) 


CSQ  has  the  capability  to  use  viscosity  which  operates  in  expansion  as 
well  as  in  compression.  This  capability  is  in  K3  and  can  be  activated  using 
variable  ARRX(12)  in  the  ARRX  array. 

In  CSQ  it  is  optional  whether  the  sound  speed  is  retained  in  the  main 
storage  arrays.  This  feature  is  retained  in  K3.  Whenever  the  sound  speed 
C$  is  not  available  for  use  in  the  viscosity  sections  of  the  program,  it  is 
approximated  by 


(2.15) 


where  C1  is  the  reference  point  and  sound  speed  for  material  i.  Here  P 
'ref 

and  p  are  cell  variables  and  5  =  5/3.  The  values  of  C-  and  5  can  be 

Vef 

changed  per  ARRX(47)  and  ARRX(48)  in  the  ARRX  array. 
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3.  Elastic-Plastic  Materials 


The  physical  relations  pertaining  to  the  elastic-perfectly  plastic  model 
in  CSQ  are  presented  in  Section  1 1-4  of  Reference  1.  Only  those  relations  that 
were  extended  for  the  K3  model  are  given  in  the  subsequent  discussion* 

The  independent  elements  of  the  stretching  tensor  are 
3  V 

dxx  =  Tf3T  (2-16) 


3V 

=  _J L 
3y 


(2.17) 


(2.18) 


1  3Vx  3Vv 

d  =  d  =  4-( — -  +  -r-i) 
xy  yx  2'3y  3x  ' 


1  3Vx  3Vz 

dxz  “  dzx  2^3z  +  3x  ^ 


l,3Vv  3Vz 

dyz  =  dzy  "  I(ir  +  W  > 


(2.19) 

(2.20) 

(2.21) 


The  volumetric  strain  rate  or  dilation  rate  is 


d  +  d 
xx  yy 


+  d 


zz 


(2.22) 
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The  stretching  deviators  are  defined  as 

A  _  .  d 
dxx  '  dxx  *  I 


dd  =  d  -  4 

yy  yy  3 


jd  _  j  d 

dZ2  -  dZZ  "  I 


dd  *  dd  =  d 
yx  xy  xy 


dd  =  dd  *  d 
zx  xz  xz 


dd  =  dd  =  d 
zy  yz  yz 


The  spin  tensor  components  are 


u  e  U  s  U  s  Q 

"xx  "yy  "zz  u 


xy 


,3V  3V 
W  s  y(  —  —  ■  ) 

yx  2v3y  3x  ' 


xz 


,  3Vy  3V_ 

W  »  ■*( — -  -  — 
zx  I'3Z  3x  ' 


i  3Vv  3V 

W  *  “  U  ■  y/  ■  .  Z  > 

yz  "zy  2vaz  3y  > 


(2.23) 

(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 


(2.32) 
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The  material  pressure  P  is 


P  =  -1(0  +0  +  a  ) 

xx  yy  zz‘ 


The  components  of  the  stress  deviator  tensor  are 


ad  =  P  -  a 
xx 


xx 


°yy  "  p  •  °yy 


4z  -  f  ■  ° 


zz 


d  _  d 
°yx  ~  axy 


xy 


Xd  a  0d 

zx  xz 


a  -  O 


xz 


0d  a  0d  a  _  0 

zy  yz  yz 


Elastic  changes  in  the  stress  deviators  are  given  by 

d?xx  _  2W  0d  -  2W  od  a  2Gdd 
*wxyxy  ^wxzxz  ^Mxx 


(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.38) 

(2.39) 


i 

l 


I 


(2.40) 
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3o 

-JM-  -  2W  ad  -  2W  ad  =  2Gdd 
yx  yx  yz  yz  yy 


(2.41) 


-  2W  crd  -  2W  od  =  2Gdd 
3t  zx  zx  zy  zy  zz 


(2.42) 


+  W  (ad  -  ad  )  -  W  od  -  W  ad  =  2Gdd 
3t  xyv  xx  yy'  yz  xz  xz  yz  xy 


(2.43) 


t  u  („d  -  „d  ) 

3t  XZV  XX  zz7 


W  0d  -  W  ad  =  2Gdd 
zy  xy  xy  zy  ^  xz 


(2.44) 


3o 


3^  +  Wyz(cTw  *  °zz]  " 


yz'  yy 


W  cxd  -  W  ad  *  2Gdd 
zx  yx  yx  zx  yz 


(2.45) 


Here  G  Is  the  shear  modulus.  Equations  2.40  -  2.45  can  be  derived  from 
equation  (6),  P.445,  Reference  2. 

The  preceding  relations  are  valid  whenever  the  material  is  in  an  elastic 
state.  Whenever  the  material  is  no  longer  in  an  elastic  state,  stress  deviators 
are  calculated  in  K3  according  to  the  preceding  relations  and  then  tested  against 
a  yielding  criteria  for  plastic  flow.  The  von  Mises  yield  criteria  used  is 


*  <dJy)2  *  <°zz>2  *  d 

*  "l/zz  *  (oxy)2  *  (<4>2 


d  0d 
xx  yy 


+  ad  od 
XX  zz 


(2.46) 
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tutor*--. 
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where  is  the  second  invariant  of  the  stress  sensor  and  Y  is  the  yield  stress. 
Whenever  this  Inequality  is  not  satisfied,  the  stress  deviators  are  adjusted  so 
that  the  stress  state  lies  on  the  yield  surface. 


Equation  (2.5)  contains  the  deviator  stress  power  P^.  It  is  given  by 


d  dd  +  ad  dd  +  ad  dd 
xx  xx  yy  yy  zz  zz 


+  2ad  dd 
xyxy 


+  2 od  dd  +  2ad  dd 
xz  xz  y z  yz. 


(2.47) 


4.  Material  Fracture 

The  K3  computation  pertaining  to  fracture  is  such  that  when  the  pressure 
P  is  less  than  a  fracture  pressure  Pfract  <  0,  void  is  introduced  into  the 
material.  As  in  CSQ,  there  are  two  K3  options  available  for  specifying  pfract- 
In  one  option  Pfract  is  a  user  defined  material  constant.  Fracture  is  based 
on  principal  stress  values  in  the  other  option. 

In  CSQ  the  maximum  and  minimum  normal  principal  stresses  and  a 2  are 

°1  ■  ?<°xx  +  V  *  *  (2’48) 

°2  ■  ?<°xx  +  V  -  ‘  «■«) 


where 


t  ■ 


(2.50) 


By  placing  a  user  specified  bound  <?fract  on  <?2  and  usin9  stress 
deviators,  Pfract  in  CSQ  is 


pfract  *  afract  +  ?  (axx  *  ayy]  +  ^  --j-**)  ♦  (ojy) 


(2.51) 


The  computation  for  P^ract  in  K3  is 


Pfract  °fract  +  Ma*  ^fxy’  ^X2*  fyz^ 


(2.52) 


where 


f  *  -1  (od  +  od  )  + 
xy  2  '  xx  yy' 


o  -  od  2  .  2 

(-^hr^  +io%) 


(2.53) 


f  =  1  (od  +  od  )  +  J(-xiLlq«)2  +  (cJd  )! 

xz  2  '  xx  zz'  " '  2  '  '  xz' 


(2.54) 


f  =  i  (od  +  ad  )  +  JuOLl^Zl)2  +  (od  \2 

yz  I  '  yy  zz'  " '  2  '  v  yz' 


(2.55) 


Note  that  equation  (2.52)  is  identical  to  equation  (2.51)  whenever 
fxy  >  *xz  and  fxy>  *yz* 
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This  fracture  model  may  be  improved  later.  Additional  discussion 
pertaining  to  fracture  is  contained  in  Section  I 1-8  of  Reference  1. 

5.  Initial  and  Boundary  Conditions 

The  K3  activity  space  is  a  region  whose  left  and  right,  bottom  and  top, 
back  and  front  sides  are  parallel  to  the  yz,  xz,  xy  coordinate  planes,  re¬ 
spectively.  The  finite  difference  mesh,  to  be  discussed  in  a  later  section, 
is  imposed  upon  this  region.  All  properties  of  the  material  within  this 
region  must  be  specified  as  initial  conditions  to  the  computational  procedure. 
Many  of  these  conditions  are  specified  via  the  equation  of  state.  The  required 
initial  input  parameters  for  a  material  are 


p  -  density 


T  or  P  -  temperature  or  pressure 


V.j  -  velocities 

o’j* j  -  stress  deviators 


along  with  some  geometrical  parameters  that  describe  the  shape  of  the  material. 
Usually  the  initial  conditions  on  the  stress  deviators  are  zero. 

Material  is  inserted  into  the  activity  space  by  defining  a  closed  shape 
around  the  material.  Several  analytic  forms  corresponding  to  various  shapes 
have  been  put  into  K3  to  aid  the  user  when  Inserting  material  into  the  mesh. 

The  geometrical  parameters  in  these  forms  must  be  specified  whenever  they  are 
used.  Figure  II-l  shows  the  five  different  shapes  currently  available  In  K3, 
along  with  the  associated  geometric  parameters  that  are  required. 


-15- 


1 .  Special  box  (edges  parallel  to  x,y,  z  axes) 


(*i i  yi  /  Z|) 

k/  ' 


N _ 

/ 

/ 


(x2.  y2/  z2) 


Required  Parameters 
x,/y,/  zi#  xJ#  y2,  z2 


2.  Tetrahedron  y>;  2]) 


x)/  yt  /  zi/  x2'  y2»  z2 

x3» y3r  Z3/  X4»  y4/  z4 


__  \  -  y4,  zd) 


(x  t  /  /\ »  zt ) 


(x2/  y2/z2) 


3.  Ellipsoid 


/  'b\ 
L-£  L- 

Tc/i  i 


(x£,  y«^  zc)  -center 
a,  b,  c  semi  major  axes  along 


x,y, z  axes 


FIGURE  11-1 


GEOMETRICAL  SHAPES  AVAILABLE  FOR  MATERIAL  REPRESENTATION 
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4.  General  solid  with  six  faces 


Required  Parameters 

»lj|»  zl,  x2r  /2r  Z2 
x3»  Xi'  z3  X4'  Y*  '  Z4 

x  s  >  y5  /  Z5f  X6  /  y**  z6 

X7'  Xr »  Z7'  x8»  ys  '  Z8 


5.  Conical  Frustrum 


(x 


1  ' 


(x2f  y2  7  z2) 


xi  /  y, ,  ^i»  x2»  y2,  Zj 

r,  r2 -radii 


FIGURE  II— 1  (continued) 


GEOMETRICAL  SHAPES  AVAILABLE  FOR  MATERIAL  REPRESENTATION 


New  material  regions  can  be  inserted  into  the  activity  space  by  combina¬ 
tions  of  insertions  and  deletions  of  the  shapes  currently  in  K3.  For  example, 
insertion  of  a  large  cube  followed  by  deletion  of  an  ellipsoid  contained 
within  the  cube  will  result  in  a  new  material  region  that  may  be  useful. 

This  capability  permits  the  user  to  insert  some  new  shapes  easily. 

The  properties  of  a  material  must  also  be  specified.  The  values  of  p,  T, 

V  ,  V  and  V  are  always  required.  Specification  of  the  stress  deviators  is 
x  y  z 

required  whenever  elastic  effects  are  to  be  treated.  Whenever  the  initial 
state  is  uniform  over  a  given  region  in  the  material,  only  a  single  value  of 
each  parameter  need  be  specified.  Fits  to  these  parameters  may  have  to  be 
supplied  whenever  these  parameters  are  dependent  on  the  x,  y  and  2  coordinates. 

K3  has  the  capability  to  impose  independent  boundary  conditions  on  the 
boundaries.  Currently  the  following  four  options  are  available  for  each 
boundary: 

0  Reflective  boundary 

0  Transmittive  boundary  -  material  can  leave  but  not  enter  mesh 

0  Transmittive  boundary  -  material  may  leave  or  enter  mesh 

0  Transmittive  boundary  -  material  may  leave  or  enter  mesh  except 

external  boundary  pressures  can  be  defined 
by  input 

Special  coding  is  usually  required  for  the  last  option. 


6.  Energy  and  Momentum  Accounting 

The  CSQ  relations  pertaining  to  energy  in  Section  11-12  of  Reference  1 
were  left  intact.  However,  modifications  to  the  momentum  accounting  relations 
were  required  as  follows. 
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The  components  of  momentum  in  the  i^  coordinate  direction  are 


Hi (t)  =  / 

1  •'u 

Vi(U)p(U)dU 

(2.56) 

where  U  is  the  computational  mesh  volume, 
has  been  added  for  the  K3  computation. 

Note  that  the  component  Hz(t) 

The  momentum 

can  be  expressed 

are  extensive  and 

lost  or  gained  due  to  material  crossing  the  mesh  boundaries 

as  several,  terms  at  the  outer  boundaries.  These  relations 
thus  will  be  expressed  in  compact  form. 

Let 

*  -  1 

if  vn>0 

=  0 

if  Vn<0 

(2.57) 

and  let  A  denote  surface  area.  Here  Vn  is  a  normal  velocity  component  to 
the  surface  A  considered.  Anchoring  these  conventions  to  the  computation 
mesh  considered,  the  following  functions  can  be  defined  for  l<.i<3, 
l<k<2: 

Ail  =  A1  eft 

if  i  =  1 

Abottom  if  i  -  2 

Aback 

if  i  =  3 

(2.58) 

A12  *  Ar1ght 

if  1  =  1 

a  A 

Mtop 

if  i  =  2 

i  A 

rtfront 

if  i  a  3 

(2.59) 

i 
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♦ik  - 1 

if  V^O 

or,  A1k 

=  0 

if  V^O 

on  A,k 

(2.60) 

The  above  conventions  can  be 

used  to  represent  the  momentum  accounting 

relations  compactly.  The  rate  that  the  i  direction 

momentum 

flows  out  of 

the  mesh  is. 

for  i^j ,  i7k,  jVk, 

*1  out  • 

/  vf  p*(2dA  -  J 

Ai2 

•  V*  n«ndA 

Ail 

+ 

/aJ2  . 

f.  viWA 

Akl 

+ 

/  »1Vk»tk2dA  - 

Ak2 

f  V-Vkp$kldA 

Akl 

(2.61) 

Similarly, 

»i  in  = 

/  v?P($i2  -  l)dA 

Ai2 

-  /A  4><*il  - 

Ail 

1 )  dA 

+ 

f  -  DdA 

j2 

-/A  W*n 

Ajl 

-  l)dA 

+ 

f.  ■  ’>  - 

Ak2 

/.  viV‘ki- 

Akl 

1 )  dA 

(2.62) 

B.  NUMERICAL  CONSIDERATIONS 
1 .  Finite  Difference  Mesh 


The  finite  difference  mesh  in  K3  is  defined  in  space  and  time.  The 
spatial  mesh  is  defined  by  the  user  while  the  time  mesh  is  usually  constructed 
by  the  program  to  fit  the  problem  considered  as  the  calculation  proceeds. 
Typically,  if  the  conditions  are  known  at  time  t^,  the  finite  difference 
algorithm  provides  the  conditions  at  a  slightly  later  time  t^.  The  time 
step  is  At  =  t^+1  -  t^  for  this  cycle.  The  next  cycle  would  have  another  time 
step  At  »  tA+2  -  t  +j.  This  process  is  continued  until  conditions  are  known 
for  the  time  range  of  interest.  Restrictions  on  the  computed  time  step  that 
pertain  to  algorithm  stability  are  presented  in  Section  B-3. 

The  K3  spatial  mesh  is  imposed  by  the  user  on  the  activity  space  considered. 
All  of  this  space  is  partitioned  into  a  three  dimensional  grid  by  three  sets  of 
parallel  planes.  These  mesh  planes  are  always  parallel  to  the  xy  or  xz  or  yz 
coordinate  planes.  The  solid  boxes  formed  by  the  intersections  of  these  mesh 
planes  are  termed  the  cells.  Since  K3  is  Eulerian  in  nature,  the  mesh  planes 
remain  fixed  throughout  the  entire  calculation.  They  are  temporarily  distorted 
during  the  Lagrangian  phase  of  a  given  cycle,  but  are  returned  to  their  original 
position  for  the  next  cycle. 

Material  and  energy  flow  through  the  fixed  mesh  and  perhaps  into  and  out  of 
the  mesh  at  the  outer  boundaries.  Boundary  condition  computations  determine 
whether  flow  is  permitted  across  the  outer  edges.  The  mesh  should  be  defined 
sufficiently  large  so  that  any  material  flow  of  significance  is  contained 
within  it. 

Let  Imax»  Jmax  and  Kmax  denote  the  number  of  mesh  increments  in  the  x,  y 
and  z  directions,  respectively,  as  determined  by  the  planes  that  partition  the 

activity  space.  Then  Imax  Jmax  Kmax  is  the  total  number  of  cells  in  the  mesh. 

The  total  nunber  of  grid  points  (xj,  y^,  zk)  is  Ipax  Jpax  Kpax, 


where 
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=  Im,„  +  1 »  +  1  and  K  _v  =  Kma>,  +  1.  A  simplified  schematic 

pax  max  pax  max  pax  max  r 

of  a  K3  mesh  is  shown  in  Figure  1 1-2 .  In  this  figure  =  6,  j  =  6  and 

max  max 

Kmax  =  thus  there  are  ^4  cells  and  245  grid  points. 


as 


It  is  apparent  from  Figure  I 1-2  that  each  cell  can  be  uniquely  identified 
|ijk|.  A  typical  cell  is  shown  in  Figure  1 1-3.  Note  the  center  of  the  cell 
with  coordinates  (x^  +  +  \/2>  zk  +  1/2^  and  center  points  of  the 

six  cell  faces.  The  volume  U^.k  of  the  cell  is 


Uijk  =  (xi+l  '  xi^yj+l  ‘  yj^zk+l  "  zk^  ^2,66) 

The  cells  in  the  mesh  need  not  all  be  of  the  same  size.  Variable  spatial 
zoning  in  the  coordinate  directions  can  be  used  to  obtain  variable  cell  size 
in  the  mesh.  Some  variable  computational  zoning  relations  are  presented  in 
Reference  1.  These  are  available  in  K3GEN. 

Material  zoning  and  rezoning  for  the  CSQ  finite  difference  mesh  is 
discussed  in  Section  1 1 1-3  of  Reference  1.  In  K3  material  rezoning  was  not 
considered.  However,  material  zoning  was  implemented  in  K3GEN.  Since  the 
CSQ  zoning  presentations  in  Section  1 1 1-3  of  Reference  1  can  readily  be 
extended  to  three  dimensions,  no  further  discussion  on  K3  material  zoning  is 
given  here. 


2.  Finite  Difference  Relations 


The  finite  difference  forms  in  K3  are  based  on  the  physical  relations 
required  for  three  dimensional  flow.  The  K3  algorithm  containing  the  dif¬ 
ference  forms  is  very  similar  in  structure  and  procedure  to  the  CSQ  finite 
difference  algorithm.  At  the  beginning  of  each  computational  cycle  the 
problem  conditions  are  assumed  known.  The  difference  equations  and  associated 


5 


*  max  I pax 
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/ 

FIGURE  11-2 

SCHEMATIC  OF  K3  SPATIAL  MESH 
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boundary  conditions  determine  the  conditions  at  a  slightly  later  time  in  the 
current  Lagrangian  coordinate  system.  The  spatial  mesh  becomes  distorted 
during  this  portion  of  the  computational  cycle.  The  Eulerian  part  of  the 
cycle  then  redefines  the  distorted  spatial  mesh  to  its  original  position. 

Other  computations  may  then  be  performed  to  complete  the  cycle. 

The  K3  finite  difference  algorithm  requires  three  spatial  indices  i,  j,  k 
and  one  time  index  i.  These  indices  can  be  used  to  indicate  where  quantities 
are  evaluated  relative  to  the  space  time  algorithm.  Assuming  that  the  indices 
i,  j  and  k  are  used  per  the  spatial  mesh  as  indicated  in  Section  Bl ,  the 
velocities  can  be  written 

Vx(i,  j  +  1/2,  k  +  1/2,  i  +  1/2) 

Vy(i  +  i/2,  j,  k  +  1/2,  i  +  1/2) 

Vz(i  +  1/2,  j  +  1/2,  k ,l  +  1/2) 

to  indicate  that  they  are  evaluated  at  the  left,  bottom  and  back  faces  of  cell 
|ijk  jat  time  t^  +  +  t&+i).  Similarly,  cell  centered  quantities 

such  as  pressure  P  can  be  written 

P(i  +  1/2,  j  +  1/2,  k  +  1/2,  i) 

to  indicate  they  are  evaluated  at  the  center  of  cell  jijk}  at  time  t^. 

The  three  dimensional  finite  difference  relations  for  the  Lagrangian 
portion  of  K3  are  extensive  and  will  not  be  presented  here.  Since  the  three 
dimensional  counterpart  of  the  CSQ  difference  forms  can  frequently  be  easily 
deduced,  the  user  can  consult  Section  IV  of  Reference  1  to  determine  informa¬ 
tion  about  the  K3  difference  forms. 
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The  Eulerian  section  of  K3  rezones  the  results  of  the  Lagrangian  portion 
of  each  cycle  so  that  the  mesh  overlays  its  position  at  the  start  of  the  cycle. 
This  requires  the  movement  of  mass,  energy  and  other  physical  quantities 
between  the  various  cells.  In  CSQ  a  geometrical  interpretation  coupled  with 
conservation  rules  is  used.  The  approach  is  one-diminsional  in  nature  with 
necessary  two-dimensional  corrections.  Further,  one  dimensional  rezoning  is 
used  for  each  of  the  i  and  j  strips  in  the  CSQ  mesh.  Biasing  due  to  this 
approach  is  minimized  by  interchanging  the  i-j  order  on  alternate  cycles. 

The  one-dimensional  approach  in  CSQ  was  retained  in  K3.  Extensions  and 
modifications  had  to  be  made  to  ensure  that  the  appropriate  corrections  were 
made  and  biasing  kept  miminal.  Note  that  on  each  cycle  there  are  three  one¬ 
dimensional  passes,  since  there  are  three  independent  coordinate  directions. 

The  order  of  the  i,  j  and  k  columns  in  K3  relative  to  the  time  cycles 
is  as  follows: 


Cycle 

- ►  0 

1 

2 

3 

4 

5 

•  •  • 

l  j 

i 

k 

j 

i 

k 

•  •  • 

Column  j 

- ►  k 

j 

i 

k 

j 

i 

•  •  • 

Order  ] 

i 

k 

j 

i 

k 

J 

•  •  • 

Note  that  the 

column  order  is  repeated 

every 

three  cycles 

.  It  is  believed 

that  permuting 

the  column  order 

in  K3 

as  shown 

here  will 

minimize  biasing 

to  the  one-dimensional  column  approach  used. 

The  Eulerian  finite  difference  relations  in  Section  V  of  Reference  1  are 
essentially  unchanged  for  K3,  since  they  are  one  dimensional  in  character. 
Thus,  they  are  not  presented  here.  Also,  the  steps  listed  in  subroutine  EULER 
for  the  alternating  direction  method  are  valid  for  K3  provided  "alternating" 
and  "strip"  are  replaced  by  "permuting"  and  "column",  respectively. 
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3.  Time  Step  Controls 


The  time  step  for  the  K3  finite  difference  scheme  must  be  selected  so 
that  the  entire  computation  remains  stable.  This  is  accomplished  by  selecting 
the  time  step  Atcs  so  that  for  any  active  cell  ijk 


At  <Mi n 
cs— 


i+1 


*ijk  + 


r5jv 


xi 


yi+l  -  yi 
Tlf^  ’  *ijk  + 


yj 


+  v 


yj+i 


where 


*1jk  =  eCs 


if 


AV  >0 

X— 


A  + 


V 


A2  +  e2C2 


if  AVx<0 


(2.67) 


(2.68) 


'"ijk  *  eCs  if  AV° 

(2.69) 

*  B  +  Vb2  +  e2C2  if  AVy<0 

n1jk  -  cCs  if  Av^° 
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A  "  BiCs  -  VVx 


(2.71) 


B  -  B£C$ 


B  AV 

q  y 


(2.72) 


B*CS  -  BqAVZ 


(2.73) 


e  *  Max 


w  -y 

i  +  v. 


(2.74) 


where  i  is  the  material  index  and  vB  is  Poisson's  ratio  of  the  zth  material. 

i 

Let  Atc$  be  the  smallest  value  that  results  from  application  of  the 
right  side  of  (2.67).  The  Courant  form  of  the  time  step  is 

Stcs  -  V~tcs  (2.75) 


where  0<f^_<l.  The  default  value  for  f^  in  K3  is  0.6. 

The  rate  of  increase  of  At  between  cycles  in  K3  is  controlled  as  in 
CSQ.  If  At^  is  the  At  of  the  last  cycle,  then  At^nc  is  defined  as 

4t1nc  -  V*'  <2-76> 

The  default  value  of  fg  is  1.068000433.  K3GEN  further  reduces  the  value 
of  At  for  the  first  cycle  of  a  computation. 


I 
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I 
I 

_  Whenever  the  time  step  is  to  remain  within  a  certain  range,  the  input 

£  parameters  At|nax  and  Atmi.n  are  used  to  provide  the  upper  and  lower  bounds. 

The  final  value  of  At  in  this  instance  is 

I 

“  *  “«  t4W  4W  <2-”) 

I 

Here  Atpad  is  idle  in  K3  since  the  radiation  option  in  CSQ  was  left 
^  intact. 

I 

I 

I 

r 

T 


I 

I 

t 

E 

£ 

I 

I 
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III  PROGRAMMING  CONSIDERATIONS 


\ 

l 


A.  RESTART  AND  PLOT  TAPE  FORMATS 

The  following  is  the  structure  used  by  K3  input/output  units  1,  2,  and 
3;  K3GEN  units  1  and  5;  and  K3PLT  input  units  1  through  5.  All  data  are 
unformatted  (binary).  The  general  file  layout  is  shown  in  Figure  III-l  below, 
followed  by  a  description  of  the  individual  records  in  each  record  block. 


Record  1 


Label 


Record  Block  2  I  EOS  Data  Block 


Record  Block  31  Cycle  Dumps  1  Through  n 


Record  Block  3, 


Record  Block  3 

n 


Figure  III-l.  Restart/Plot  Tape  File  Structure 


« 
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Record  1  -  Problem  Identification  Format  (10A10) 

(LABU),  4=1,  10) 

Record  Block  2  -  EOS  Data  Block 


Record  2a 

a 

(ND(i),  4=1,6) 

Record  2. 

D 

(A(t) ,  4  =  1,  ND(1 ) ) 

Record  2C 

(K(4),  4=1,  ND(2) ) 

Record  2d 

(A(4),  4  =  1,  ND(3))  ‘ 

There  are  ND (4 )  of  the  following  record  [ND(4)>_  0] 

Record  2n 
e 

(A(4),  4=1,  ND(5) ) 

Record  2f 

Following  record  is  present  only  if  ND(6)  =  1 
(A(4),  4  =  1  ,  8267) 

where  ND(1)  _<  10 

ND(2)  =  142 
ND (3 )  =  7321 
0  _<  ND(4)  <  10 
ND(5)  =  6851 
ND (6 )  =  0  or  1 

Record  Block  3  -  Cycle  Dump 

Record  3a  NI,  NF,  (NARRX(4) ,  4  =  1,  NI),  (ARRX(i) ,  4  =  1,  NF) 

where  NI  =  112 
NF  =  1168 

For  the  following  record,  IMAX  =  NARRX(13) 

JMAX  =  NARRX (14) 

KMAX  =  NARRX (104) 

Record  3b  ((VXL(4,m),  4=1,  JMAX),  m  -  1,  KMAX), 

((VYL(4,m),  ft  -  1,  IMAX),  m  -  1,  KMAX), 

( (VZL(4,m) ,  4*1,  IMAX),  m  «  1,  JMAX) 


flgwwyy  ..... .  -.*»,  .....  _  ..  ...... 
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Record  3.  There  are  KMAX*IMAX  records  of  the  form 
(STORE(a) ,  *  =  1 ,  JNVAR) 

where  JNVAR  =  NARRX(20).  Each  record  is  an  ik-column,  i.e. 
all  cell  variables  for  a  column  of  cells  1  <  j  <  JMAX. 

Following  the  last  column  record  there  is  the  start  of  the  next  cycle 
dump  (Record  3 J  or  an  end-of-file. 

a 
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B.  CELL  VARIABLES  AND  STORAGE  ARRAYS 

A  set  of  variables  must  be  retained  for  each  cell  in  the  K3  mesh.  The 
storage  variables  are  defined  in  such  a  way  that  only  those  variables  re¬ 
quired  for  the  problem  at  hand  are  retained.  Table  III-l  lists  the  variables 
that  may  be  associated  with  each  cell.  The  basic  number  of  variables  required 
is  ten  and  may  range  up  to  40.  To  facilitate  the  storage  of  a  variable  number 
of  cell  variables,  an  array  called  LK(t)  is  defined  which  contains  the  relative 
location  of  each  cell  variable  associated  with  a  particular  cell.  An  *  value 
is  associated  with  each  cell  variable,  e.g.,  t  =  1  is  pressure,  i  =  2  is  tem¬ 
perature,  etc.  LIU  and  LK(i)  are  equivalenced  in  the  coding  of  K3  for  ease 
of  computations.  Example:  If  there  are  NVAR  variables  retained  per  cell 
and  all  data  are  stored  sequentially  in  an  array  called  STORE,  it  follows  that 
the  temperature  of  the  third  cell  is  located  at  ST0RE(2*NVAR+LK2). 

Since  there  are  NVAR*IMAX*JMAX*KMAX  cell  variables,  it  would  be  impracti¬ 
cal  to  store  all  of  this  information  in  central  memory  except  for  very  small 
values  of  IMAX,  etc.  or  very  large  computers.  Thus,  cell  variables  are  stored 
in  ECS  (or  LCM)  to  the  extent  of  available  ECS,  with  excess  information  stored 
in  Random  Acess  Mass  Storage  (disk).  In  ECS,  the  cell  variables  are  arranged 
in  terms  of  a  index  n  given  by 

"tjk  -  Ww  *  k-'IJmx  +  i 

where  i,  j,  k  are  the  indices  of  a  particular  cell.  The  variables  for  a  given 
cell  i,  j,  k  are  then  located  in  ECS  locations  [ ( n-1 ) *NVAR+1 ]  to  [n*NVAR], 
assuming  the  cell  variables  are  located  in  ECS  and  not  on  disk.  If  some,  or 
all,  cell  variables  are  located  on  disk,  matters  become  more  complicated.  In 
this  case.  Random  Mass  Storage  indices  are  required  to  locate  the  desired  cell 
information  on  disk  (all  variables  associated  with  any  one  cell  comprise  one 
record  on  disk).  These  indices  are  also  stored  in  ECS.  Thus,  ECS  must  be 
partitioned  to  accommodate  as  many  cell  variables  as  possible  and  also  the 


TABLE  III-l 


Cell  Storage  Variables 


Variable 

z 

LKU)  Val 

p 

(pressure) 

1 

l 

T 

(temperature) 

Cs 

(sound  speed) 

3 

★ 

U 

(cell  volume) 

4 

* 

Vy 

(y  velocity,  bottom  of  cell) 

5 

3 

Vx 

(x  velocity,  left  of  cell) 

6 

4 

E. 

(•  specific  energy) 

7 

5 

M. 

(•  mass) 

8 

6 

Ex 

(x  specific  energy) 

9 

7 

Mx 

(x  mass) 

10 

8 

$x 

(x  material  volume  fraction) 

11 

9 

•v 

(void  volume  fraction) 

12 

* 

xr 

(mean  free  path) 

13 

★ 

0^ 

XX 

(stress  deviator) 

14 

★ 

0^ 

xy 

(stress  deviator) 

'  15 

* 

0^ 

y.y 

(stress  deviator) 

16 

* 

HE  i 

variable  1 

17 

★ 

HE  i 

variable  2 

18 

★ 

reserved  for  code  expansion 

19 

* 

position  for  start  of  material  arrays 

20 

★ 

Vz 

(z  velocity,  back  of  cell) 

21 

★  * 

0^ 

xz 

(stress  deviator) 

22 

★ 

•J. 

(stress  deviator) 

23 

★ 

°iz 

(stress  deviator) 

24 

★ 

reserved  for  code  expansion 

25 

★ 

♦Indicates  variables  retained  only  If  required  for  problem  under  consid¬ 
eration.  Actual  values  for  a  given  computation  can  be  determined  from  the  listing 
of  the  NARRX  array  output  In  K3GEN  and  K3  runs.  Note  that  LK(1)  =  NARRX(55),  etc. 

+<T2  =ST0RE(JVAR+LK20+1 ) 

E2  "STORE ( JVAR+LK20+2 ) 

M2  "STORE (JVAR+LK20+3) 

?3  "STORE (JVAR+LK20+4) 


**  V,  is  required  in  every  three 
c  dimensional  problem.  Also, 
LK21=LK20+3(N-1 )+l  where  N 
is  the  number  of  materials. 


Mn.1  "STORE (JVAR+LK20+ [N-2J ) 
JVAR  "{(cell  number  -  1)NVAR} 


indices  to  Random  Mass  Storage  records  for  the  cell  variables  that  must  be 
stored  on  disk.  The  partitioning  is  done  by  an  alternate  entry  to  Subrou¬ 
tine  GETPUT,  and  is  based  upon  the  amount  of  available  ECS,  information 
supplied  by  Subroutine  MEM.  It  is  advantageous  to  try  to  fit  all  of  the 
cell  variables  into  ECS  by  setting  the  size  of  ECS  as  large  as  possible 
and/or  scaling  down  the  number  of  cells,  since  using  disk  results  in  a 
longer  running  program  (the  ratio  of  wall  clock  time  to  central  processor 
time  increases  considerably). 

The  actual  location  of  cell  variables,  either  on  disk  or  in  ECS,  is 

transparent  to  the  program  and  user  (except  where  run  time  is  increased) 

since  the  storing  and  fetching  of  these  variables  is  handled  by  one  routine, 

GETPUT.  This  routine  works  with  one  "strip"  at  a  time  where  a  strip  is  all 

of  the  variables  associated  with  a  column  or  row  of  cells,  i.e.,  all  cells 

where  i  and  k  are  constant  and  j  varies  from  1  to  j  ,  or  i  and  j  are  con- 

max 

Stant  and  k  varies  from  1  to  k  .  etc.  To  handle  the  cell  variables 

max 

associated  with  a  strip  a  series  of  STORE  arrays  (STORE,  STOREL,  STORER,  etc), 
dimensioned  NVAR*  max  ( imax »  Jmax.  kmax),  are  defined  in  K3GEN  and  K3*  GETPUT 
will  fetch  a  strip  from  ECS  or  disk  to  one  of  the  STORE  arrays,  the  main  rou¬ 
tines  will  process  and  modify  information  where  required,  and  then  call  GETPUT 
to  store  the  information  back  to  ECS  or  disk  before  moving  on  to  the  next  strip. 

Additional  arrays  are  required  to  hold  the  cell  location  and  size  informa¬ 
tion  (X,  Y,  Z,  AX,  aY,  aZ)  and  also  the  velocities  at  the  top,  right,  and  front 
of  the  mesh  (VYL,  VXl,  VZL).  These  latter  arrays  are  also  stored  in  ECS.  At 

present  the  dimensions  of  the  arrays  are  set  to  accommodate  up  to  80  cells  in 

each  X,  Y  or  Z  direction.  If  these  are  increased,  variable  NECSL  in  K3GEN 

must  be  redefined  to  indicate  the  new  length  of  common  /VEL/,  and  KMAXM, 

JMAXM,  IMAXM,  also  in  K3GEN,  must  be  reset  to  the  maximum  number  of  cells 
al 1  owed . 

Two  arrays,  NARRX  and  ARRX,  are  located  at  the  beginning  of  blank 
common  In  K3GEN  and  K3.  These  are  user-accessible  arr^s  with  provisions 
made  In  K3  for  changing  any  variable  In  either  array  by  input  cards  In  K3. 

Table  1 1 1-2  lists  the  variables  and  their  use  (-—  indicates  an  unused 
variable)  in  these  two  arrays. 
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TABLE  II 1-2 

Blank  Common  Variable  List  for  K3 


NARRX 

Other  Names 

Purpose,  Quantity,  Section 
or  Equation  of  Definition 

1 

Number  of  Variables  in  ARRX  Array  Before  Geometry 
Arrays  (114  for  present  code) 

2 

First  Location  in  ECS/LCM  for  Main  Mesh  Storage 
Arrays 

3 

— 

4 

— 

5 

Flag  for  Suppression  of  Thermal  &  Plasma  Conducts 
0  -  Normal 

1  -  Suppress 

6 

HE  Burn  Flag  (IBURNHE  input  variable) 

0  -  No  Burn 

+  -  Burn  HE  (Idle  -  not  operational) 

7 

Current  Number  of  Energy  Flow  Subcycles  Allowed 

8 

1 

Maximum  Number  of  Energy  Flow  Subcycles  Allowed 

9 

Energy  Flow  Cycle  Number 

10 

j 

1 

11 

I 

Mesh  Index 

12 

J 

i 

Mesh  Index 

13 

IMAX 

1 

See  Section  II-B-1 

14 

JMAX 

See  Section  II-B-1 

15 

I  PAX 

See  Section  II-B-1 

16 

JPAX 

See  Section  II-B-1 

17 

I  ACT 

See  Section  II-B-1 

18 

JACT 

See  Section  II-B-1 
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NARRX 

Other  Names 

Purpose,  Quantity,  Section 
or  Equation  of  Definition 

19 

NVAR 

Number  of  cell  variables 

20 

JNVAR 

JMAX*NVAR 

21 

NORAD 

Energy  Flow  Switch 

0  -  No  Energy  Flow  Computation 

1  -  Energy  Flow  Computation  on  (Idle  -  not  opert 

22 

NVOID 

Cell  Void  Flag 

0  -  Retain  Cell  Void  Fractions 

1  -  Do  Not  Retain 

23 

NEXTV1 

Print  Variable  Flag 

24 

NEXTV2 

Print  Variable  Flag 

25 

NEXTV3 

Print  Variable  Flag 

26 

ICYCLE 

Hydro  Cycle  Number 

27 

I COUNTS 

Short  Edit  Number 

28 

ICOUNTL 

Long  Edit  Number 

29 

NDUMP 

Restart  Dump  Number 

30 

MDUMP 

Plot  Dump  Number 

31 

ICYCLES  ! 

End  of  Problem  Cycle 

32 

NEXTP 

Flag  for  Special  Edit 

33 

NDTCS 

Cell  Number  Controlling  Atcs 

34 

NOTCSN 

Same  as  33  for  Next  Cycle 

35 

NDTRAD 

1 

1 

Cell  Number  Controlling  Atratj  (Idle) 

36 

i 

i 

NOTRADN 

Same  as  35  for  Next  Cycle  (Idle) 

37 

IGM 

Geometry  Flag  6 

0  =  Rectangular 

1  =  Cylindrical  (Idle  -  not  operational) 
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Other  Names 


Purpose,  Quantity,  Section 
or  Equation  of  Definition 


NUMMT 

NUMMX 


Number  of  Material 

Number  of  X  Materials  (NUMMT-1) 


KAB(1 ) 


KAB(l)  to  KAB(IO)  are  Hollerith  Fields  with  Problem 
Identification 


KAB(2) 

KAB(3) 

KAB(4) 

KAB(5) 

KAB(6) 

KAB(7) 

KAB(8) 

KAB(9) 

KAB(IO) 

KLAB(l) 


Current  Date 
Current  Clock  Time 

KLAB(l)  to  KLAB(4)  are  Hollerith  Fields  with 
Edit  Titles 


KLAB ( 2 ) 
KLAB(3) 
KLAB (4) 
KPLOT 


Short  Edit  Line  Printer  Plot  Flag 
0  -  No  Plot 
1  -  Material  Plot  Only 

1  -  Momentum  Plot  Only 

2  -  For  Both  Plots  (Normal) 
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NARRX 

Other  Names 

Purpose,  Quantity,  Section 
or  Equation  of  Definition 

55 

LKl.LK(l) 

LK1  to  LK25,  See  Table  III-l 

56 

LK2,LK(2) 

57 

LK3,LK(3) 

58 

LK4,LK(4) 

59 

LK5,LK(5) 

60 

LK6,LK(6) 

61 

LK7,LK(7) 

62 

LK8,LK(8) 

63 

LK9,LK(9) 

64 

LK10,LK(10) 

65 

LK11 ,LK(11 ) 

66 

LK1 2  >LK(  12) 

67 

LK13,LK(1 3) 

68 

LK14,LK(14) 

69 

LK1 5,LK( 15) 

70 

LK16,LK(16) 

71 

LK17,LK(17) 

72 

LK18,LK(18) 

73 

LK19,LK(19) 

74 

LK20,LK(20) 

75 

LK21 ,LK(21 ) 

76 

LK22,LK(22) 

- 

77 

LK23,LK(23) 

78 

LK24,LK(24) 

79 

LK25,LK(25) 

NARRX 


Other  Names 


Purpose,  Quantity,  Section 
or  Equation  of  Definition 


— 

80 

NELEP 

Elastic  Plastic  Switch 

0  -  No  EP  Calculation 

1  -  EP  Calculation  on 

81 

NPORFG(l) 

Distended  Material  Switch  for  Material  1 
=  0  not  distended 
=  1  can  be  distended 

82 

NP0RFG(2) 

\ 

83 

NP0RFG(3) 

84 

NP0RFG(4) 

1 

85 

NP0RFG(5) 

f 

86 

NP0RFG(6) 

\Same  as  NPORFG(l)  except  for  material  indicated 

87 

NP0RFG(7) 

( 

88 

NP0RFG(8) 

89 

NP0RFG(9) 

90 

NPORFG(IO) 

i 

/ 

Tracer  Particle  Storage  Index  (Idle  -  not  operational) 

91 

MMMM(l) 

92 

MMMM(2) 

Tracer  Particle  Storage  Index  (Idle  -  not  operational) 

93 

MMMM(3) 

Maximum  I  Mesh  Size  Possible  (Size  of  Dimension 
Statements,  80  for  Present  Code) 

94 

MMMM(4) 

Same  as  93  Except  for  J  Mesh 

95 

MMMM(5) 

I  Mesh  Index  for  Activation  Calculation 

96 

MMMM(6) 

J  Mesh  Index  for  Activation  Calculation 

97 

MMMM(7) 

Current  Number  of  Energy  Flow  Subcycles 

98 

MMMM(8) 

— 

99 

MMMM(9) 

— 

100 

MMMM(IO) 
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ARRX 


Other  Names 


Purpose,  Quantity,  Section 
or  Equation  of  Definition 


1 


Minimum  Velocity  for  Cell  Activation 


2 


Minimum  Fractional  Temperature  Difference  for  Cell 
Activation 


3 

4 


Temperature  Convergence  Factor  for  Energy  Balance 
Iterations 


5 


6 

7 


8 

9 

10 


11 

12 


13 

14 

15 


Flag  to  Suppress  Writing  of  Global  Quantities  on 
History  Tape  55.  =  0  Write  Global  Quantities 
=  1  Do  Not  Write 

Temperature  Lower  Limit  for  EOS  Iterations 

Flag  for  Viscosity  Form.  If >0  Use  Q  in  Compression 
Only.  «  0  Use  Normal  Q  in  Both  Expansion  and 
Compression.  <  0  Use  Normal  Q  in  Compression  and 
Linear  Term  in  Expansion 

Temperature  Iteration  Limit  Factor 


16 


Other  Names 


Purpose,  Quantity,  Section 
or  Equation  of  Definition 


17 

... 

18 

— 

19 

— 

TEMPBG 

Temperature  Assigned  to  Void  Cell  if  NORAD 

21 

— 

22 

— 

23 

— 

24 

— 

25 

— 

26 

Factor  in  Energy  Flow  Time  Step  Control 

27 

Factor  in  Opacity  Average 

28 

Factor  in  Energy  Flow  Time  Step  Control 

29 

Factor  in  Opacity  Average 

30 

TMASS 

Total  Mass  in  Mesh 

31 

TMASSD 

Total  Dot  Mass  in  Mesh 

32 

TMASSX 

Total  X  Mass  in  Mesh 

33 

TIME 

t  =  Problem  Time 

34 

DT 

At 

35 

DTH 

At' 

36 

DTOLD 

At  of  Last  Cycle 

37 

DTCS 

Atcs 

38 

DTCSN 

Atcs  for  Next  Cycle 

39 

DTRAD 

Atrad 

ARRX 

Other  Names 

Purpose,  Quantity,  Section 
or  Equation  of  Definition 

40 

DTRADN 

Atrad  for  Next  Cycle 

41 

DTMAX 

At 

max 

42 

DTMIN 

Atmin 

43 

FRACDT 

fl 

44 

DTINCR 

f2 

45 

BQ 

Bq 

46 

BL 

B£ 

47 

QCS 

Ciref 

48 

QGAM 

5 

49 

VCUT 

Cut  Off  Momentum  for  Short  Edit  Momentum  Plots 

50 

ETOT 

initial  +  (see  E<?s-  (2.106)(2. 107)  in  Ref. 

51 

HUP 

V*> 

52 

HRAO 

nx(t) 

53 

HFOR 

Hz(t) 

54 

ERADIN 

Radiation  Energy  Into  Mesh 

55 

ERADOUT 

Radiation  Energy  Out  of  Mesh 

56 

EINTIN 

Internal  Energy  Into  Mesh 

57 

E INTOUT 

Internal  Energy  Out  of  Mesh 

58 

EKIN 

Kinetic  Energy  Into  Mesh 

59 

EKOUT 

Kinetic  Energy  Out  of  Mesh 

60 

HUP  IN 

Hy(t)  Into  Mesh 

61 

HRAO IN 

Hx (t )  Into  Mesh 

Hy(t)  Out  of  Mesh 


62 


HUPOUT 


ARRX 


Other  Names 


Purpose,  Quantity,  Section 
or  Equation  of  Definition 


63 

HRADOUT 

64 

HFORIN 

65 

HFOROUT 

66 

AMDIN 

67 

AMDOUT 

68 

AMXIN 

69 

AMXOUT 

70 

BSWTOP 

71 

BSWSID 

72 

BSWBOT 

73 

BSWINS 

74 

BSVIFOR 

75 

BSWBAC 

76 

RADTOP 

77 

RADS ID 

78 

RADBOT 

79 

RADINS 

80 

RADFOR 

81 

RADBAC 

82 

TSTOP 

83 

GMINES 

84 

GMAXES 

85 

BURNP1 

Hx(t)  Out  of  Mesh 
Hz(t)  Into  Mesh 
Hz(t)  Out  of  Mesh 
Dot  Mass  Flow  Into  Mesh 
Dot  Mass  Flow  Out  of  Mesh 
X  Mass  Flow  Into  Mesh 
X  Mass  Flow  Out  of  Mesh 

Hydro  Boundary  Condition  Option  (Input  Parameters) 

Same  as  70 
Same  as  70 
Same  as  70 
Same  as  70 
Same  as  70 

Energy  Flow  Boundary  Condition  Option  (Input  Parameter) 

Same  as  76 

Same  as  76 

Same  as  76 

Same  as  76 

Same  as  76 

Current  Problem  Stop  Time 

Ymin 

Ymax 

D/N  -  HE  Burn  Parameter 
HE  Predetonation  Pressure 


86 


BURNP2 


Other  Names 


Purpose,  Quantity,  Section 
or  Equation  of  Definition 


ARRX 


87 

TSQUEN 

88 

PFR114 

89 

BWORK 

RADK1 

91 

RADK2 

92 

RADK3 

93 

RADK4 

94 

RADK5 

95 

RADK6 

96 

RADK7 

97 

RADK8 

98 

RADK9 

99 

PIE 

100 

TWOPIE 

101 

AMINM 

102 

DVOLVM 

103 

EVOLCK 

104 

VCHECK 

105 

VOLCK 

106 

RADAVS 

107 

EMPTY 

-  I 

Total  Energy  from  Sources 

Parameter  Used  for  Extrapolation  of  Boundary  Pressures 

4  as/c 
4  o$/(3c) 

16  a$/(3c) 

16  as/c 

as/6 

as/16 

as/12 

as 
2/3 

tt  (  =  3.1415926..) 

2n 

Round  Off  Test  for  Empty  Cell 

Cut  Off  Test  Parameter  for  Kinetic  Energy  Corrector 

Maximum  Relative  Charge  in  Volume  Allowed  in 
Energy  Balance 

Smallest  Allowed  Velocity 

Round  Off  Test  Parameter  for  Volume  Checks 

Minimum  Fractional  Difference  Between  Adjacent  Zone 
Temperatures  to  Compute  Boundary  Mean  Free  Path  Values. 

If  0  -  Always  Use  Cell  Center  Values 

Scale  Factor  for  Kinetic  Energy  Corrector  (Normally  1) 
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ARRX 

Other  Names 

Purpose,  Quantity,  Section 
or  Equation  of  Definition 

108 

EARRX(l) 

109 

EARRX(2) 

110 

EARRX(3) 

Extra  Array  Not  Used  in  Current  Coding 

111 

EARRX(4) 

112 

EARRX(5) 

113 

EARRX(6) 

114 

EARRX(7) 

•  A.  CONTROL  CARDS  FOR  K3  PROGRAMS 

The  run  sequence  for  a  K3  calculation  begins  with  PREK3.  This  program 
generates  the  COMDECKs  necessary  for  compilation  of  K3GEN,  K3,  and  K3PLT. 

The  procedure  used  in  the  following  examples  of  control  card  setup  was  to 
punch  the  COMDECKs,  then  physically  insert  the  punched  cards  in  the  proper 
place  in  the  K3GEN,  etc.,  job  decks.  This  procedure  can  certainly  be  modi¬ 
fied  to  suit  the  needs  or  desires  of  individual  users,  limited  only  by  the 
capabilities  of  the  UPDATE  program. 

1 .  PREK3 

Job  card 
ATTACH,  A,  K3. 

UPDATE  (P  =  A,  C  =  PREK3,  Q) 

RETURN,  A. 

FTN  (I  =  PREK3,  L  =  0) 

RETURN,  PREK3 
LGO. 

REWIND,  DKCSQ. 

COPYCF  (DKCSQ,  PUNCH)  Punch  deck  for  K3GEN,  K3,  and  K3PLT. 

REWIND,  DKCSQ.  IOUT  input  variable  set  to  1  to 

COPYCF  (DKCSQ,  PUNCH)  direct  output  to  file  DKCSQ. 

REWIND,  DKCSQ. 

COPYCF  (DKCSQ,  PUNCH) 

EXIT. 

7/8/9 
*C  PRECSQ 
7/8/9 

-  data  card  for  PREK3  - 

7/8/9 

6/7 /8/9 


f 
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2.  K36EN 

Job  card 

ATTACH  (OLDPL ,  K3) 

UPDATE  (P,  Q) 

RETURN,  OLDPL. 

FTN  (I,  LCM  =  I,  L  =  0) 
RETURN,  COMPILE. 

REQUEST,  TAPE1  ,  *PF. 

LGO  (PL  =  50000) 

CATALOG  (TAPE1 ,  RESTART) 
EXIT. 

7/8/9 

*AF  INPUT,  YANK$$$ 

-  COMDECKs  created  by  PREK3 
*C  CSQGEN 

7/8/9 

-  data  cards  for  K3GEN  - 
7/8/9 

6/7 /8/9 


3.  K3 

Job  card 

ATTACH  (OLDPL,  K3) 
UPDATE  (P,  Q) 

RETURN,  OLDPL. 

FTN  (I,  LCM  -  I,  L  ■  0) 
RETURN,  COMPILE. 

ATTACH  (TAPE1 ,  RESTART) 
STAGE  (TAPE 2,  POST) 
LDSET  (PRESET  =  ZERO) 
LGO  (PL  -  50000) 

EXIT. 
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3.  K3  (continued) 

7/8/9 

*AF  INPUT,  YANK$$$ 

-  COMDECKs  created  by  PREK3  - 
*C  CSQ,  EP 

7/8/9 

-  data  cards  for  K3  - 
7/8/9 

6/7 /8/9 

Note  the  compilation  of  the  EP  deck  which  contains  the  elastic-plastic 
routines.  This  deck  may  or  may  not  be  required  depending  upon  the  problem 
being  run.  Other  decks  which  may  be  compiled  such  as  the  History  package. 
Tracer  package,  etc.,  are  not  operational  in  the  main  K3  code. 

4.  K3PLT 

Job  card 

ATTACH  (OLDPL ,  K3) 

UPDATE  (P,  Q,  R) 

REWIND,  COMPILE. 

FTN  (LCM  =  I ,  I,  L  =  0) 

UPDATE  (F,  P,  C  =  SCORS,  R) 

REWIND,  SCORS. 

FTN  (I  *  SCORS,  L  =  0,  B  =  XLIB,  LCM  =  I ) 

RETURN,  SCORS,  COMPILE,  OLDPL. 

REWIND,  XLIB. 

LIBEDT. 

RETURN,  XLIB. 

ATTACH  (ARCLIB,  ARCLIB,  ID  =  PRDLIB) 

LIBRARY  (SLIB,  ARCLIB) 

LDSET  (PRESET  »  ZERO) 

LGO. 

EXIT. 


0 

D 
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4.  K3PLT  (continued) 

7/8/9 

*AF  INPUT,  YANK$$$ 

-  COMDECKs  from  PREK3 

*C  CSQPLT 

7/8/9 

7/8/9 

LIBRARY  (SLIB,  NEW  =  4000) 

ADD  (*,  XL  IB) 

FINISH 

ENDRUN 

7/8/9 

7/8/9 

6/7/8/9 

The  modified  SCORS  package  of  routines  is  put  into  a  library  in  order 
to  reduce  the  size  of  the  executable  program  since  only  those  routines  re¬ 
quired  are  loaded.  No  data  cards  are  required  since  inputs  are  defined 
interactively,  i.e.,  through  keyboard  inputs  as  the  program  executes. 
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B.  INPUT  INSTRUCTIONS 


1 .  PREK3 

The  program  PREK3  is  a  preprocessor  for  K36EN,  K3  and  K3PLT.  It 
generates  COMDECKs  for  use  by  the  UPDATE  processor.  It  must  be  used  to 
produce  a  program  acceptable  to  the  FTN  compiler.  Only  one  card  is  read 
by  PREK3.  Its  format  is  as  follows: 


Card  1 

Format  (415) 

d-5) 

UMAX  *  largest  mesh  dimension  (IMAX,  JMAX  or  KM  AX) 
20  <  UMAX  _<  80 

(6-10) 

NVAR  =  maximum  number  of  variables  per  cell 

10  _<  NVAR  <  40 

(11-15) 

MAXMS  =  maximum  number  of  materials 

2  _<  MAXMS  <  10 

(6-20) 

I OUT  *  output  switch 

=  0,  punched  card  output 
=  1,  output  to  file  DKSCQ 

PREK3  also  produces  for  user  edification  a  listing  of  expected  program 
memory  requirements  based  upon  the  generated  COMMON  sizes.  This  listing  has 
not  been  updated  from  CSQ  requirements  (K3  requires  more  space),  hence  should 
not  be  used  to  size  programs  until  such  time  as  it  is  modified. 

2.  K3GEN 

The  program  K3GEN  generates  a  new  problem  or  rezones  a  previous  K3 
calculation.  The  following  suggestions  are  included  to  aid  the  user.  Their 
meaning  will  be  clearer  after  going  through  the  input  Instructions. 

1.  EOS  and  high  explosive  data  should  be  checked  with  the  test 
program  CKE0S2. 

2.  ANEOS  options  -1,  0,  and  1  should  not  be  used  with  K3. 

3.  The  cell  activator  works  in  Increasing  coordinate  direction. 

Initial  action  should  be  in  the  lower  left  back  corner  of  the 
mesh  for  calculational  speed  considerations. 


4.  During  material  zoning,  subcell  centers  are  tested  to  determine 
if  a  subcel  1  is  inside  a  zoning  region. 

5.  If  more  than  one  material  zoning  region  incloses  a  given  spatial 
region,  the  first  region  data  is  used. 

6.  During  a  rezone  any  card  inputs  are  processed  before  the  old  mesh 
data  are  considered.  Because  of  5  above,  this  means  the  card  in¬ 
puts  will  override  old  mesh  data. 

7.  Cell  centers  are  tested  to  determine  if  a  cell  is  inside  a  region 
defined  as  high  explosive. 

8.  ECS/LCM  required  is  discussed  in  the  text  of  report. 

9.  Only  ANEOS  options  3  and  4  and  tabular  EOS  information  should  be 
used  for  porous  materials. 

10.  During  a  rezone  the  ECS/LCM  required  is  the  larger  of  that  required 
for  the  old  or  new  mesh. 

11.  Rules  fo'*  automatic  determination  of  number  of  cell  variables  (NVAR) 
using  data  on  card  4. 

Start  with  NVAR  *  10 

Add  1  if  NORAD  >  0 

Add  1  if  NVOID  =  0 

Add  1  if  CSSAVE  *  1 
Add  1  if  VOLSAVE  =  1 

Add  6  if  NELEP  =  1 

Add  2  if  IBURNHE  t  0* 

Add  3  for  each  material  over  two. 


Tape  Units  K3GEN 

1.  Restart  output  tape  (for  K3) 

2.  Rezone  Input  tape  (from  K3)* 

5.  Input  tape  from  CHDCSQ"^  or  internal  storage  during  rezones 

12.  EOS  input  tape  for  tabular  data 
72.  Internal  Storage 

_ 77.  Random  Mass  Storage 

*  Not  operational  In  present  K3  code, 
t  Not  tested 
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a.  K3GEN  Inputs 


Card  1.  Format  (8A10) 

80  Column  Problem  Identification 
Any  BCD  Information 


Card  2.  Format  (1615) 

Variable  1.  NUMMT 


Variable  2.  NREZONE 


Variable  3.  NCHDZN 


Variables  NEOS(KEOS) 

[KEOS  »  1,..., NUMMT] 


Number  of  materials 

Cannot  be  larger  than  defined  in  PREK3 
2  _<  NUMMT  _<  10. 

=  0  for  generation  of  new  problem 

>  0  for  rezone  of  previous  calculation 
with  old  EOS  data  using  the  NREZONE— 
dump  on  input  unit  2.* 

<  0  for  rezone  of  previous  calculation 
with  new  EOS  data  (define  as  for  new 
problem)  using  the  j NREZONE |^dump  on 
input  unit  2.  If  NREZONE  t  0,  insert 
Card  Set  R  after  Card  Set  3  and  see 
discussion  at  the  end  of  this  section.* 

=  0  No  input  data  from  CHDCSQ 
=  1  To  read  input  tape  from  CHDCSQ  on 
tape  Unit  5.’'  If  NCHDZN  *  1,  then  code 
sets  NREZONE  =  0  above;  i.e.,  both 
features  cannot  be  used  at  same  time. 

Equation  of  state  number  for  material  KEOS 

>  0  tabular  (on  Tape  12) 

*  -1,  -2,..., -NUMMT  (ANEOS) 


*  Not  operational  In  present  K3  code, 
t  Not  tested 


Any  inputs  for  the  ANEOS  package  go  here.  The  formats  are  given  in  Appendix 
F  of  SLA-73-0477  (Reference  4)  or  section  b,  following  this  section.  This  is 
the  same  as  Card  Set  12  for  CHART  D.  The  identification  variable  in  the  ANEOS 
input  format  should  agree  with  NEOS(KEOS)  above.  There  is  one  set  of  cards 
for  each  unique  NEOS(KEOS)  <  0.  They  can  be  in  any  order. 

Omit  this  card  set  if  NREZONE  >  0  on  Card  1. 


Card  Set  R. 

If  NREZONE  /  0  (Card  1)  insert  Card  Set  R  here.  See  discussion  at  end  of 
K3  instructions. 


Card  4. 

Variable  1 . 
0-5) 


Variable  2. 
(6-10) 


Variable  3. 
(11-15) 


Variable  4. 
(16-20) 


Format  (1615) 
IGM 


NORAD 


NVAR 


NVOID 


-  geometry  switch 
=  0  rectangular 

=  1  cylindrical  (not  operational) 

-  energy  flow  switch 
=  0  no  energy  flow 

>  0  energy  flow  with  NORAD  the  maximum 
number  of  subcycles  allowed.  (Not 
operational ) 

=  number  of  variables  stored  per  cell. 

=  0  for  code  to  calculate  (normal) 

<  0  for  code  to  calculate  and  then 
add  -NVAR  to  result 

-  void  switch 

=  0  code  carries  cell  void  fractions 
(normal).  Use  0  if  any  solids  are 
present  In  mesh. 

*  1  code  does  not  carry  void  fractions. 
Use  1  only  if  voids  can  never  be  present. 
(\very  high  temperature,  gases,  etc.) 


Variable  5 
(21-25) 


CSSAVE 


-  switch  to  save  sound  speed  for  printing 
=  0  do  not  save  (normal) 

=  1  save  sound  speed 


Variable  6. 
(26-30) 


Variable  7. 
(31-35) 


Variable  8. 
(36-40) 


Variable  9. 
(41-45) 


Variable  10. 
(46-50) 


Variable  1 1 . 
(51-55) 


Variable  12. 
(56-60) 


♦See  Table  III 
°xx’  °xy’  °yy 


V0LSAVE 


NEXT  VI 


NEXTV2 

NEXTV3 

fCLEP 

NPORQUS 

KATMOS 


-  switch  to  save  cell  volume  for  printing 
=  0  do  not  save  ( normal ) 

=  1  save  volume 

-  print  switch  for  nonstandard  edit 
variable  in  3rd  from  last  edit  column, 
printed  variable  is  ST0RE(LK(NEXTV1 ) ) .* 
Default  values  are: 

11  (4  )  if  NORAD  =  0 

*  x/ 

13  Ur)  if  NORAD  f  0 

-  same  as  NEXTV1  except  for  second  from 
last  edit  column.  Default  value  is  7(E.). 

-  same  as  NEXTV1  except  for  last  edit 
column.  Default  value  is  9(Ex). 

-  switch  for  elastic  plastic  calculation 
=  0  EP  calculation  off 

*  1  EP  calculation  on 

-  switch  for  porous  material  calculation 
=  0  porous  material  calculation  off 

=  1  porous  materials  present- input  card  7. 

not  used 


1.  Note  if  NEXTV1 ,  2,  3  ■  14, 

instead  of  o^  ,  . 

xx’  xy’  yy 


15,  or  16  printed  variable  is 
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If  ,■ 


Variable  13.  IBURNHE 
(61-65) 


Variable  14.  NTRACE 
(66-70) 


Card  5.1 

Variable  1 . 
(1-10) 


Variable  2. 
(11-20) 

Variable  3. 
(21-30) 


Format  (8E1Q.3) 


=  0  if  no  high  explosives  are  to  be 
burned 

=  1  if  high  explosives  are  present 
(not  operational). 

=  -1  same  as  +1  except  with  predetona¬ 
tion  print  in  CSQ  (tells  when  cell 
starts  burning). 

If  IBURNEHE  i  0,  see  Card  Set  11. 

(Not  operational ) 

-  switch  for  tracer  particles 
=  0  no  tracer  particles 
=  1  input  tracer  particle  data  on 
Card  Set  12.  (Not  operational) 

=  -1  retain  tracer  positions  from 
rezone  data.  (Not  operational) 

=  2  retain  tracer  positions  from 
rezone  data  and  input  more  particle 
data  on  Card  Set  12.  (Not  operational) 

LAST  TWO  OPTIONS  APPLY  ONLY  TO  REZONES 
(NREZONE  t  0  ON  CARD  2). 


-  X  position  for  active  mesh  if 
XACT  <  XI  NS  IDE  (see  below), 
code  sets  XACT  *  • 

-  Y  position  for  active  mesh  if 
YACT  <  YBOTTOM,  code  sets  YACT  *  - 

-  Z  position  for  active  mesh 
If  ZACT  j<  ZBACK 

code  sets  ZACT  *  - 


Variable  4. 

XI  NS  I DE 

-  X(l) 

(31-40) 

Variable  5. 

YBOTTCM 

-  Y(D 

(41-50) 

Variable  6. 

ZBACK 

=  Z(l) 

(51-60) 

Card  5.2 

Format  (8E10.3) 

Variable  1 . 
(1-10) 

BL 

-  linear  viscosity  coefficient 

Variable  2. 
(11-20) 

BQ 

-  quadratic  viscosity  coefficient 

If  BL  +  BQ  =  0,  Code  Sets  BL  =  0.1 
and  BQ  =  2.0. 

Variable  3. 
(21-30) 

TIME 

-  initial  problem  time  (normally  0). 

If  a  rezone  is  being  done  (NREZONE  t  0 
on  Card  2),  input  -  1 .  to  retain  time 

from  rezone  data. 

Variable  4. 
(31-40) 

TEMPBG 

-  background  temperature  (used  as 
initial  temperature  in  regions  not 
included  in  Card  Set  10.)  Default 
value  is  0.02567785ev  (298°K). 

Card  6. 

Format  (8F5.0. 

E10.3) 

Variable  1 . 

BSWTOP 

=  mass  boundary  switch  for  top  of  mesh 

(1-5) 

(options  given  below). 

Variable  2. 
(6-10) 

BSWBOT 

=  mass  boundary  switch  for  bottom  of  mesh 

Variable  3. 
(11-15) 

BSWSID 

=  mass  boundary  switch  for  right  of  mesh. 

Variable  4 
(16-20) 


BSWINS 


=  mass  boundary  switch  for  left  of  mesh. 
(Not  used  if  IGM  =  1) . 


Variable  5. 
(21-25) 

BSW  B  AC 

=  mass  boundary  switch  for  back  of  mesh. 

Variable  6. 
(26-30) 

BSW  FOR 

=  mass  boundary  switch  for  front  of  mesh. 

Variable  7. 
(31-35) 

RAOTOP 

=  radiation  boundary  switch  for  top  of  mesh. 

Variable  8. 
(36-40) 

RADBOT 

=  radiation  boundary  switch  for  bottom  of 

mesh. 

Variable  9. 
(41-45) 

RADS  ID 

=  radiation  boundary  switch  for  right  of 

mesh. 

Variable  10. 

RADINS 

=  radiation  boundary  switch  for  left  of 

(46-50) 

mesh.  (Not  used  if  IGM  =  1). 

Variable  11. 
(51-55) 

RADBAC 

*  radiation  switch  for  back  of  mesh. 

Variable  12. 
(56-60) 

RADFOR 

=  radiation  switch  for  front  of  mesh. 

Variable  13. 
(61-65) 

DTSTART 

=  fraction  of  normal  time  step  to  be 
used  on  cycle  1.  If  DTSTART  _<  0,  code 

sets  DTSTART  *  0.02. 

BSW  —  Options 

RAD 

—  Options  (Not  Operational) 

0.  Reflecting  boundary 

0. 

Reflecting  boundary 

1.  Transmittive  boundary  no  mass 

1. 

Transmittive  boundary  with  mass  outside 

al  1  owed  to 

enter  mesh. 

of  mesh.  Properties  of  external  mass  are 

2.  Same  as  1 

except  mass  can  be 

interpolated  from  the  properties  Inside  mesh. 

drawn  into 

the  mesh. 

2. 

Same  as  Type  1  except  external  temperature 
is  defined  by  input  cards  In  K3.  (See  Card 
Set  14  in  CSQ  input.  Reference  1). 

3.  Same  as  2  except  that  exter 
nal  boundary  pressure  can 
be  applied.  Will  require 
special  coding  in  K3  sub¬ 
routine  PBTRP. 


3.  Transmittive  boundary  without  mass  out 
side  of  mesh.  External  temperature  is 
defined  by  input  cards  in  K3  (See  Card 
Set  14  in  CSQ  input.  Reference  1). 


Card  7.  Format  (1015) 

This  set  of  cards  is  present  only  if  NP0R0US  =  1  on  card  4.  There  is  one 
input  parameter  for  each  material.  1  =  1,  NUMMT. 

Variable  NPORFG(KEOS)  =  switch  for  porous  material 

[KE0S=1  ,...,NIMMT]  cal cul ation  for  material  KEOS. 

=  0  not  porous 

=  1  material  can  be  distended 


Zoning  the  Problem 

Each  problem  is  zoned  in  two  parts.  There  is  a  spatial  card  set  (X,  Y,  Z) 
which  defines  the  space  fixed  mesh.  A  material  definition  card  set  then  follows 
to  place  the  material  in  this  mesh.  Each  card  carries  an  identification  field 


1  to  5  to  aid 

Identification 

in  ordering  the  punch  cards. 

Concerns 

40' s 

X  spatial 

50' s 

Y  spatial 

60' s 

Z  spatial 

70'  s 

Insert  material  into  mesh. 

80's 

Delete  material  (in  relation  to  70  cards). 

90's 

Material  properties. 

The  numbers  40,  50,  60,  70  and  90  are  only  used  as  stop  cards.  Numbers  41,  93, 
etc.,  are  employed  to  define  information.  As  will  be  explained  below,  the 
order  of  the  cards  must  be 
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41  -  can  be  more  than  one. 

40 

51  -  can  be  more  than  one. 

50 

61  -  can  be  more  than  one. 

60 

71  to  74 
79 

81  to  84  -  optional 

89  -  optional 
91  to  94  -  optional 

90 
70 

NOTE:  If  present.  Card  Set  G  follows  the  90  card.  See  discussion  at  end 
of  K3GEN  instructions. 


Repeat  for  each  material  package 


Card  Set  8.  X  Zoning 

The  zoning  starts  at  the  left  of  the  mesh  and  works  to  the  right.  Note 
that  XINSIDE  =  X ( 1 )  is  defined  on  Card  5.1. 


Card  8.1  Format  (215,  4E10.3) 

Variable  1.  Integer  41 

(1-5) 


Variable  2. 
(6-10) 

Variable  3. 
(11-20) 


N  =  number  of  zones  to  be  included  with  this 
card  set. 

DX  =  width  of  first  zone  (left  side)  in  this  region 


Variable  4.  X  =  total  width  of  this  region. 

(21-30)  1 


Variable  5. 
(31-40) 


R  *  ratio  of  adjacent  zone  widths. 
(  *  AXi+1/aXi ) 
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Variable  6.  DXL  =  width  of  last  zone  (right  side)  in  this  region. 

(41-50) 


Note:  Define  the  variables  by  one  of  the  following  four  options. 
Leave  others  blank. 


Option 

N 

DX 

xt 

R 

DXL 

1 

/ 

✓ 

/ 

2 

/ 

✓ 

/ 

3 

✓ 

/ 

/ 

4 

/ 

(N>2) 

/ 

/ 

/ 

R  is  not  constant  for  Option  4 


Keep  repeating  these  cards  until  the  entire  X  mesh  is  defined.  When  finished 
insert  Card  8.2. 

Card  8.2  Format  (15)  X  zoning  stop  card 

Variable  1.  Integer  40 

(1-5) 

Card  Set  9.1  Y  Zoning 

The  zoning  starts  at  the  bottom  of  the  mesh  and  works  upward.  Note  that 
YBOTTOM  =  Y ( 1 )  is  defined  on  Card  5.1. 

The  information  on  the  two  forms  of  input  cards  is  the  same  as  Card 
Set  8  except  that  51  and  50  are  the  identification  numbers. 

Card  Set  9.2  Z  Zoning 

The  zoning  starts  at  the  back  of  the  mesh  and  works  forward.  Note  that 
ZBACK  3  Z(l)  Is  defined  on  Card  5.1. 

The  information  on  the  two  forms  of  input  cards  is  the  same  as  Card 
Set  8  except  that  61  and  60  are  the  identification  numbers. 


Material  Zoning 


In  CSQ  the  material  is  placed  in  the  space  mesh  by  defining  geometrical 
shapes  in  the  X,  Y  plane  and  the  properties  of  the  material  in  this  region. 
Shapes  can  be  inserted  and  then  partially  deleted  to  obtain  oddly  configured 
material.  For  example  the  shape 


This  method  of  material  zoning  for  CSQ  was  generalized  to  three  dimen¬ 
sions  for  K3.  The  shapes  in  K3  that  can  be  inserted  and  then  partially 
deleted  to  obtain  irregular  configured  material  are  as  follows. 

1.  Rectangular  parallelepiped  or  box  (edges  parallel  to  x,  y,  z  axis) 


Code  No .  t 
71  -  insert 


79/89  Card  Variable 


4 


X 

Y 


2 

2 


81  -  delete 


5 
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Code  No. 

72  -  insert 
82  -  delete 


3.  Ellipsoid 


Code  No.  79/89  Card  Variables 

73  -  insert  1  X0  4a 

83  -  delete  2  YQ  5b 

3  Zq  6  c 
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4.  Irregular  solid  figure  bounded  by  quadrilaterals 


(*$  -yt>h) 


Code  No. 

74  -  insert 
84  -  delete 


12 


24 


5.  Frustum  of  cone 


Code 

No. 

79/89  Card 

Variables 

75  - 

insert 

1 

X1 

5 

V2 

85  - 

del ete 

2 

Y1 

6 

Z2 

3 

Z1 

7 

rl 

4 

X2 

8 

r2 

Note:  This  will  require  two  79/89  cards. 

Material  Properties  Fits 

In  most  cases,  the  material  properties  within  a  given  geometrical  shape 
will  be  uniform  in  the  initial  zoning.  Special  cases  might  be  encountered 
which  require  variation  within  the  shape.  Provisions  are  available  to  treat 
this  situation  though  special  coding  will  normally  be  required. 

If  this  special  fit  is  used  (NFIT  =  1  on  Card  10.1)  insert  cards  below 
Instead  of  Cards  10.7  to  10.10.  Nonna! ly  used  with  shape  option  No.  III. 


T  -  C3  ♦  C45 
Vy  -  C5  (Y  -  Y0) 
Vx  -  C6  (X  -  x0) 
Vz  *  C7  (2  -  Z„) 


Variable  Following  Card 
Identification  Number 


Constant 


1,  2,  3 


C5’  C6’  C7 


No  stress  deviators  are  defined  (94  card),  i.e.. 


3  0^  =  3  =o^  =o^  3  0,  with  Fit  No.  1 

zz  yz  xz  xx  xy  yy 


Card  Set  10 

Card  10.1 

Variable  1 . 
(1-5) 


Format  (515,  A55) 
Shape  Code  (71  to  74) 


Insertion  of  material  into  mesh 
Insert  Material  Option 


Variable  2. 
(6-10) 


=  Number  of  sti>divisions  in  cell  co¬ 
ordinates  for  zoning.  (1  to  100  de¬ 
pending  on  resolution  required).  There 

3 

are  NSUB  subcells  per  cell.  The  center 
of  each  sJbcell  is  tested  against  the 
zoning  package. 


I 
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=  Number  of  delete  package  (80  maximuii) 
to  follow  and  applied  against  this 
insert  package. 

=  Material  type  indicator 

0  =  void  0  would  be  used  to  define  a 
hot  spatial  region,  i.e., 
hot  void 

>0  for  material  KEOS.  (Material  KEOS  uses 
EOS  number  NEOS(KEOS)  as  on  card  2.) 

=  0  for  uniform  region  material  properties. 
_>  1  =  Fit  number  for  nonuniform  material 
properties  (see  previous  material  proper¬ 
ties  fits  disdussion). 

Two  fits  are  presently  coded  as  examples 

(NFIT  *  1  and  707) 

for  comment. 

Card  10.2  Format  (15,  5X,  6E10.3)  Shape  Dimensions 

Variable  1.  Integer  79 

(1-5) 

Variable  2  to  7.  Shape  variables  1  to  6  (See  previous  Shape  Options  Section) 
(11-70) 

Card  10.3  Format  (15,  5X,  2E1Q.3)  Continuation  of  last  card 

Used  only  with  shape  option  IV. 

Variable  1.  Integer  79 

0-5) 

Variables  2  and  3.  Shape  variable  7  and  8 
(11-30) 


Variable  3.  NDEL 

(11-15) 


Variable  4.  KEOS 

(16-20) 


Variable  5.  NFIT 

(21-25) 


Variable  6.  Any  BCD  information 

(26-80) 
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Note:  Cards  10.4  -  10.6  are  present  only  if  NDEL  >  0  on  Card  10.1.  If 
NOEL  >  0,  there  are  NDEL  sets  of  Card  10.4  to  10.6. 

Card  10.4  Format  (15)  Delete  for  last  insert  package 

Variable  1.  Shape  Code  (81  to  84) 

(1-5) 

Cards  10.5  and  10.6 

Same  as  Card  10.2  and  10.3  except  code  is  89  instead  of  79. 


Note:  See  Material  Properties  Fits  for  Cards  10.7  to  10.10  if  NFIT  on  Card 
10.1  is  greater  than  zero.  Any  of  the  cards  with  identification  num¬ 
ber  91  to  94  can  be  left  out  if  a  zero  or  default  input  is  desired. 
They  can  be  in  any  order.  Card  10.11  (identification  number  90)  must 
always  be  present. 


Card  10.7 

Format  (15.  5X.  7E10.3) 

Density 

Variable  1 . 
(1-5) 

Integer  91 

Variable  2. 
(11-20) 

p  =  material  density 

Card  10.8 

Format  (15.  5X.  7E10.3) 

Temperature  or  Pressure 

Variable  1 . 
(1-5) 

Integer  92 

Variable  2. 

T  =  material  temperature  [if  T<0 
(298°K)] 

OR 

code  sets  T  =  0.02567785  ev 

Variable  1 . 
(1-5) 

Integer  92 

Variable  2. 
(11-20) 

o 

P  *  material  pressure  (dynes/cm  ) 

Variable  3. 
(21-30) 

1.  (flag) 
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Card  10.9 

Format  (15.  5X.  7E10.3) 

Variable  1 . 

Integer  93 

(1-5) 

Variable  2. 

V  =  X  velocity  of  material 

(11-20) 

X 

Variable  3. 

V  =  Y  velocity  of  material 

(21-30) 

y 

Variable  4. 

V,  =  Z  velocity  of  material 

(31-40) 

z  J 

Card  10.10* 

Format  (15.  5X.  7E10.3) 

Variable  1 . 

Integer  94 

(1-5) 

Variable  2. 

c* 

(11-20) 

XX 

Variable  3. 

°d„ 

(21-30) 

yy 

Variable  4. 

d 

a,. 

(31-40) 

zz 

Variable  5. 

d 

a  „ 

(41-50) 

xy 

Variable  6. 

<4 

(51-60) 

yz 

Variable  7. 

cd 

(61-70) 

xz 

♦Card  10.10 

should  only  be  used  if  NELEP  =  1  or 

Card  10.11 

Format  (15) 

Variable  1 . 

Integer  90 

0-5) 

X,  V  and  Z  Velocities 


Stress  Deviators 


Stop  card  (must  be  present) 


If  NF IT  «  707  (Card  10.1),  insert  Card  Set  G  here.  See  discussion  at  end 
of  K3GEN  instructions. 


To  define  more  material  regions  go  to  Card  10.1  and  start  over.  When  finished 
Insert  Card  10.12. 


-73- 


Card  10.12 

Variable  1 . 
(1-5) 


Format  (15) 
Integer  70 


Stop  card 


Card  Set  G 


Stable  configuration  with  gravity 


This  set  of  cards  is  used  to  define  a  stable  configuration  in  a  gravi¬ 
tational  field  (Y  direction  only).  The  information  is  used  by  the  material 
properties  calculation  for  NFIT  =  707  (Card  10.1  in  K3GEN).  When  present, 
this  card  set  follows  the  90  card  (Card  10.11  in  K3GEN).  There  can  be 
several  sets  of  these  cards;  one  for  each  material  insert  package  desired. 


If  there  are  several  sets  of  these  cards,  care  must  be  used  to 

top  and  bottom  pressures. 

Card  G.l  Format  (15,  5X,  6E10.3) 

Variable 

0-5) 

1. 

Integer  707 

Variable 

(11-20) 

2. 

YB  =  smallest  Y  for  package 

Variable 

(21-30) 

3. 

PB  =  pressure  at  YB+ 

Variable 

(31-40) 

4. 

YT  =  largest  Y  for  package 

Variable 

(41-50) 

5. 

PT  =  pressure  at  YT+ 

Variable 

(51-60) 

6. 

T  For  a  constant  temperature  atmosphere 

temperature.  If  T  =0,  code  set  T  =  , 
For  an  altitude  dependent  temperature 
negative  number  and  Card  Set  G. 2. 

Variable 

(61-70) 

7. 

RHOG  =  approximate  density  at  boundary  witl 
pressure  for  starting  iterations. 

t  Oefine  either  PB  or  PT  and  leave  other  blank 
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Card  Set  G. 2  Format  (15,  5X,  2E1Q.3)  Temperature  -  Altitude  Relation 

Present  only  if  T  <0  (Variable  6  on  Card  G.l). 

Insert  one  card  for  each  desired  point.  Start  at  bottom  of  mesh  and  work 
upward.  Stop  with  Card  G. 3.  V 

Variable  1.  Integer  701 

d-5) 

Variable  2.  Y  For  following  temperature. 

(11-20) 

Variable  3.  Temperature 

(21-30) 


Card  G.3  Format  (15)  Stop  Card 

Present  only  if  Card  Set  G.2  is  used. 

Variable  1.  Integer  700 

(1-5) 


The  equation-of-state  input  data  information  presented  here  is  repro¬ 
duced  from  Reference  1  for  the  convenience  of  the  user.  This  information 
can  also  be  found  in  Appendix  F  of  Reference  4. 

The  card  set  forms  the  necessary  data  cards  for  Card  Set  3  in  K3GEN, 
Card  Set  12  in  CHART  13,  Card  Set  2  in  CKE0S2  and  any  other  program  which 
uses  the  ANEOS  equation-of-state  package.  Note  that  there  are  8  optional 
forms  with  4  different  input  formats.  The  first  variable  on  the  first  card 
of  each  form  must  agree  with  the  data  on  Card  2  of  K3GEN. 

Inputs  for  ANEUS  options  0  to  +4 

Card  1.  Format  (13,  15.  12.  5A10.  2E10.3) 

Equation-of-state  number  (negative  number). 

-1  to  -10  for  K3GEN 

Library  equation-of-state  number  if  desired; 
otherwise  zero.+ 

Used  only  with  a  library  equation-of-state.  This 
variable  determines  the  type  of  analytic  calculation 
(see  variable  2,  card  2  below).  If  out  of  range  0  to 
4,  or  library  information  is  only  for  a  gas,  this  input 
is  ignored. 

Fifty-column  identification  label  :  any  BCD  information 

RHUG  -  The  initial  density  for  the  Hugoniot  calculation 
If  zero,  the  calculation  is  skipped.  If  negative,  the 
initial  density  is  taken  to  be  the  reference  density 
(variable  3,  card  2  below). 


Variable  1 . 

(1-3) 


Variable  2. 
(4-8) 


Variable  3. 
(9-10) 


Variable  4-8. 
(11-60) 


Variable  9. 
(61-70) 


See  Appendix  C  in  R3  for  contents. 
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Variable  10. 
(71-80) 


Cards  2,  3,  and 
In  the  listing. 


Variable  1 . 
(1-10) 


Variable  2. 
(11-20) 


Variable  3. 
(21-30) 


Variable  4. 
(31-40) 


Variable  5. 
4 1 -SO) 


•  •'  it)  f  6. 


THUG  -  The  initial  temperature  for  the  Hugoniot 
calculation.  If  zero,  the  calculation  is  skipped. 

If  negative,  the  initial  temperature  is  taken  to  be  the 
reference  temperature  (variable  4,  card  2  below). 

If  a  library  equation  of  state  is  requested, 
no  further  data  cards  are  required. 

4^  Format  (8E10.3) 

the  following  variables  are  called  ZB(I),  I  =  1,  24. 

The  number  of  elements  in  this  material. 


Switch  for  type  of  equation-of-state.  See  Reference  8 
for  these  options. 

0.  -  Solid-gas  without  electronic  terms  and  without 
detailed  treatment  of  the  liquid-vapor  region. 

1.  -  Solid-gas  with  electronic  terms  but  without 

detailed  treatment  of  the  liquid-vapor  region. 

2.  -  Gas  only  with  electronic  terms. 

3.  -  Same  as  0.,  but  with  a  detailed  treatment  of 

the  liquid-vapor  region. 

4.  -  Same  as  1.,  but  with  a  detailed  treatment  of  the 

liquid-vapor  region. 

pQ  -  Reference  density. 

Tq  -  Reference  temperat^e 

If  T  _<  0,  code  sets  Tq  =  0.02567785  ev  (298°K). 

PQ  -  Reference  pressure  (normally  0). 

B0  -  Reference  bulk  modulus  (positive  number) 
or 

(-SQ)  -  Constant  in  linear  Hugoniot  shock-particle 
velocity  relation  (negative  number) 


i 

i 
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Variable  7.  r  -  Reference  Gr'uneisen  coefficient. 

(61-70)  0 


Variable  8. 
(71-80) 


Variable  9. 
(1-10) 


Variable  10. 
(11-20) 


Variable  11 . 
(21-30) 


±0q  -  Reference  Debye  temperature. 

If  6  =0,  code  sets  6  =  0.025. 
o  o 

If  0Q  >  0,  use  high  temperature  approximation 
If  0q  <  0,  calculate  complete  Debye  functions 

-  Parameter  for  selection  of  model . 

T r  =  - 1 ,  SI ater  theory, 

=  0,  Dugdale  and  MacDonald  theory, 

Tp  =  1 ,  free-volume  theory, 
or 

-  Constant  in  linear  Hugoniot  shock -particle  velocity 
relation.  Input  variable  is  defined  in  relation 

to  variable  6. 

30^  -  Three  times  the  limiting  value  of  the  Gr’uneisen 
coefficient  for  large  compressions,  usually  either 
2  or  0. 

When  a  value  of  2  is  used,  =  2/3. 

E$  -  Zero  temperature  separation  energy. 


Variable  12. 
(31-40) 


Variable  13. 
(41-50) 


Tm  -  Melting  temperature  (positive  number) 
or 

(-Em)  -  Energy  to  the  melting  point  at  zero  pressure 
from  the  reference  point  (negative  number). 

-  Parameter  for  low  density  Pc  modification  to 
move  critical  point  (normally  zero). 
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Variable  14. 
(51-60) 


Variable  15. 
(61-70) 


Variable  16. 
(71-80) 


Variable  17. 
(1-10) 


Variable  18. 
(11-20) 


Variable  19. 
(21-30) 


Variable  20. 
(31-40) 


Variable  21 . 
(41-50) 


Variable  22. 
(51-60) 


C54  -  Parameter  for  low  density  Pc  modification  to 
move  critical  point  (normally  zero). 

If  =  0  and  ^  0,  codes  sets  Cg4  =  0.95. 

H  -  Thermal  conductivity  coefficient.  If  zero,  thermal 
o  J 

conduction  is  not  included.  Note  that  the  units  of 

r 

H  =  H  T  41  are  ergs/ (cm  sec  eV). 

-  Temperature  dependence  of  thermal  conduction 
ccoefficient  (see  variable  15). 

pmin  “  Lowest  a^owed  solid  density,  usually  about  0.8pQ. 
If  zero  or  negative,  code  sets  pmin  =  0.8  pQ. 

Parameter  D-j 

Parameter  Dg  I 

Parameter  0^  V 

Parameter  \ 

Parameter  Dc  i 
5  / 


Solid-solid  phase  transition 
parameters  (normally  0). 


Variable  23.  -  Heat  of  fusion  to  determine  melt  transition  parameters. 

(61-70)  If  =  0,  no  transition  is  included. 

If  <  0,  code  sets  Hf  =  1.117  x  1012  T^/Afergs/gm) , 
where  A  is  the  average  atomic  weight. 

Note :  Code  will  run  slower  if  the  melt  transition  is 
included.  Use  only  when  necessary  and  after 
testing. 


Variable  24 


p4/ps  -  Ratio  of  liquid  to  solid  density  at  melt  point 
or 

(-p  )  -  Density  of  liquid  at  melt  point, 
or 

1  +  aV  =  1  +  V^  -  Vg  -  Change  in  volume  at  melt  plus  one. 
Note:  In  the  first  option,  the  input  number  is  between  0 
and  1;  in  the  second  it  is  negative;  and  in  the  third 
greater  than  one. 

If  i  0  and  pa/ps  =  0,  code  sets  pa/p$  =  0.95. 

For  a  gaseous  equation  of  state  (type  2),  variables  5  to  14  and  17  to  24 
are  read  but  not  used. 


Card  5  Format  (5(F5.0,  E10.3)) 

There  is  one  set  of  the  following  variables  for  each  element  in  variable  1, 
card  2. 

1=1,  number  of  elements. 

Variable  Odd.  Z(I)  -  Atomic  number  of  element. 

Variable  Even.  Unnormalized  atomic  number  fraction  of  element 

[COT (I )]. 
or 

-  (Unnormalized  atomic  weight  fraction  of  element.) 

All  elements  should  be  defined  in  the  same  way. 

Inputs  for  ANE05  Option  -1 

Card  1.  Format  (13.  15.  12.  5A10.  2E10.3) 

Variable  1.  Equation-of-state  number  (negative  nunber) . 

(1-3) 


Variable  2 
(4-8) 


Library  equation-of-state  number  if  desired; 
otherwise  zero. 


Variable  3. 
(9-10) 


Variable  4-8 
(11-60) 


Variable  9 
(61-70) 


Variable  10. 


Cards  2,  3.  and 

In  the  listing. 

Variable  1 . 
(1-10) 


Variable  2. 
(11-20) 


Variable  3. 
(21-30) 


Variable  4. 
(31-40) 


Used  only  with  a  library  equation-of-state.  This 
variable  determines  the  type  of  analytic  calculation 
(see  variable  2,  card  2  below). 

Must  be  -1  in  this  case. 

Fifty-column  identification  label  :  any  BCD  information. 

RHUG  =  The  initial  density  for  the  Hugoniot  calculation. 
If  zero,  the  calculation  is  skipped.  If  negative,  the 
initial  density  is  taken  to  be  the  reference  density 
(variable  3,  card  2  below). 

THUG  =  The  initial  temperature  for  the  Hugoniot  calcu¬ 
lation.  If  zero,  the  calculation  is  skipped.  If  nega¬ 
tive,  the  initial  temperature  is  taken  to  be  the  ref¬ 
erence  temperature  (variable  4,  card  2  below). 

If  a  library  equation-of-state  is  requested, 
no  further  data  cards  are  required. 

_4  Format  (8E10.3) 

the  following  variables  are  called  ZB(I),  I  =  1,  24. 

The  number  of  elements  in  this  material. 

Switch  for  type  of  equation  of  state. 

Mist  be  -1  for  this  option. 

pQ  -  Reference  density. 

Tq  -  Reference  temperature. 

If  T0  £  0»  code  sets  TQ  =  0.02567785  eV  (298°K). 


Variable  5 
(41-50) 


Blank 
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Variable 

(51-60) 


Variable 

(61-70) 


Variable 

(71-80) 


Variable 

d-10) 


Variable 

(11-20) 


Variable 

(21-30) 


Variable 

(31-40) 


Variable 

(41-50) 


Variable 

(51-60) 

Variable 

(61-70) 


6.  SQ  -  Reference  point  bulk  sound  speed. 

7.  r  -  Reference  Gruneissen  coefficient. 

8.  Blank. 

9.  Sj  -  Constant  in  linear  Hugoniot  shock-particle 
velocity  relation. 

or 

(-100)  -  For  power  series  representation  of  Pu 

n 

Enter  variables  18  to  21 
or 

(-S-|  -100)  -  For  power  series  representation  with 
Sj  for  low  density. 

10.  Blank. 

11.  C  -  Heat  capacity. 

If  1  0,  code  sets  Cy  =  3  No.  k. 

12.  Tm  -  Melting  temperature  (positive  number) 

or 

(-Em)  -  Energy  to  the  melting  point  at  zero 
pressure  from  the  reference  point  (negative 
number) . 

13.  Blank. 

14.  Blank. 

15.  Hq  -  Thermal  conductivity  coefficient.  If  zero, 

thermal  conduction  is  not  included.  Note  that 
the  units  of  H  *  H0T^41  are  ergs/(cm  sec  eV). 
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Variable  16. 
(71-80) 


Variable  17. 
(1-10) 


Variable  18. 
(11-20) 


Variable  19. 
(21-30) 


Variable  20. 
(31-40) 


Variable  21 . 
(41-50) 


Variable  22. 
(51-60) 


Variable  23. 
(61-70) 


Variable  24. 
(71-80) 


Card  5 

There  is  one  set 
card  2. 

1  »  1 ,  number  of 
Variable  Odd. 
Variable  Even. 


C^i  -  Temperature  dependence  of  thermal  conduction 
coefficient  (see  variable  15). 

pmin  "  '-owest  allowed  solid  density,  usually  about 
0.8pQ.  If  zero  or  negative,  code  sets  pmin  =  0.8pQ. 


Constants  in  power  series  representation 
of  PH. 

Enter  only  if  variable  9  ^  -100. 

Blank  otherwise. 


Blank. 

Blank. 

Blank. 


Format  (5(F5.0,  E10.3) ) 

of  the  following  variables  for  each  element  in  variable  1, 

elements. 

Z(I)  -  Atomic  number  of  element. 

Unnormalized  atomic  number  fraction  of  element 
[COT (I )]. 
or 

-(Unnormalized  atomic  weight  fraction  of  element.) 
All  elements  should  be  defined  in  the  same  way. 


Variable  1 . 
0-3) 


Equation-of-state  number  (negative  number). 


Variable  2. 
(4-8) 

Variable  3. 
(9-10) 


Variable  4-8. 
(11-60) 

Variable  9. 
(61-70) 


Variable  10. 
(71-80) 


Cards  2,  3,  and 

In  the  listing. 

Variable  1 . 
(1-10) 

Variable  2. 
(11-20) 

Variable  3. 
(21-30) 


Library  equation-of-state  number  if  desired;  otherwise  zero. 

Used  only  with  a  library  equation-of-state.  This  variable 
determines  the  type  of  analytic  calculation  (see  variable  2, 
card  2  below).  Must  be  -2  in  this  case. 

Fifty-column  identification  label:  any  BCD  information. 

RHUG  =  The  initial  density  for  the  Hugoniot  calculation. 

If  zero,  the  calculation  is  skipped.  If  negative,  the 
initial  density  is  taken  to  be  the  reference  density 
(variable  3,  card  2  below). 

THUG  =  The  initial  temperature  for  the  Hugoniot  calculation. 

If  zero,  the  calculation  is  skipped.  If  negative,  the  initial 
temperature  is  taken  to  be  the  reference  temperature 
(variable  4,  card  2  below). 

If  a  library  equation-of-state  is  requested, 
no  further  data  cards  are  required. 

_4  Format  (8E10.3) 

the  following  variables  are  called  ZB( I ) ,  I  «  1,  24. 

The  number  of  elements  in  this  material. 

Switch  for  type  of  equation-of-state. 

Must  be  -2  for  this  option. 

pQ  -  Reference  density. 


T0  -  Reference  temperature. 

If  T0  i  °.  code  sets  T0  =  0.02567785  eV  (298°K). 

K,  -  Constant  in  Rosseland  opacity  expression. 

O 

Units  are  cm  /gm. 

Blank. 

Variable  7.  (Y  -  1)  -  Specific  heat  ratio  minus  one. 

(61-70) 

K2  -  Constant  in  Rosseland  opacity  expression  Units  are  eV. 
If  <2  <  0,  code  sets  «2  =  .025  eV. 

Blank. 

Variable  10.  Blank. 

(11-20) 

Cy  -  Heat  capacity. 

If  Cv  <  0,  code  sets  Cy  =  3/2  No  k. 

Blank. 


Variable  13.  Blank 

(41-50) 

Variable  14.  Blank. 

(51-60) 

Variable  15.  Blank. 

(61-70) 

Variable  16.  Blank. 

(71-80) 

Variable  17.  Blank. 

(1-10) 


Variable  1 1 . 
(21-30) 


Variable  12. 
(31-40) 


Variable  8. 
(71-80) 


Variable  9. 
(1-10) 


Variable  4. 
(31-40) 


Variable  5. 
(41-50) 


Variable  6. 
(51-60) 


I 


I 


1 

Variable  18. 

(11-20) 

Blank 

I 

Variable  19. 

(21-30) 

Blank 

I 

Variable  20. 

(31-40) 

Blank 

I 

Variable  21 . 
(41-50) 

Blank 

I 

Variable  22. 

(51-60) 

Blank 

I 

Variable  23. 

(61-70) 

Blank 

I 

Variable  24. 

(71-80) 

Blank 

Card  5. 


Format  (5(F5.0.  E10.3)) 


r 

i 

im 


>  T 


[ 

c 


I 

I 


There  is  one  set  of  the  following  variables  for  each  element  in  variable  1, 
card  2. 

1=1,  number  of  elements. 


Variable  Odd. 
Variable  Even. 


Z(I)  -  Atomic  number  of  element. 

Unnormalized  atomic  number  fraction  of  element  [COT ( I )D » 
or 

-(Unnormal ized  atomic  weight  fraction  of  element.) 

All  elements  should  be  defined  in  the  same  way. 

Inputs  for  ANEOS  Option  -3 

Card  1.  Format  (13.  15.  12.  5A10,  2E10.3) 

Variable  1.  Equation-of- state  number  (negative  number). 
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Variable  2. 
(4-8) 

Variable  3. 
(9-10) 


Variables  4-8 
(11-60) 

Variable  9. 
(61-70) 


Variable  10. 
(71-80) 


Cards  2,  3.  and 

In  the  listing. 

Variable  1 . 
(1-10) 

Variable  2. 
(11-20) 

Variable  3. 
(21-30) 

Variable  4. 
(31-40) 


Library  equation-of-state  number  if  desired; 
otherwise  zero. 

Used  only  with  a  library  equation-of-state.  This 
variable  determines  the  type  of  analytic  calcula¬ 
tion  (see  variable  2,  card  2  below).  Must  be  -3 
in  this  case. 

Fifty-column  identification  label:  any  BCD  information. 

RHUG  -  The  initial  density  for  the  Hugoniot  calculation. 
If  zero,  the  calculation  is  skipped.  If  negative,  the 
initial  density  is  taken  to  be  the  reference  density, 
(variable  3,  card  2  below). 

THUG  =  The  initial  temperature  for  the  Hugoniot  calcu¬ 
lation.  If  zero,  the  calculation  is  skipped.  If 
negative,  the  initial  temperature  is  taken  to  be  the 
reference  temperature  (variable  4,  card  2  below). 

If  a  library  equation-of-state  is  requested, 
no  further  data  cards  are  required. 

1  Format  (8E10.3) 

the  following  variables  are  called  ZB(I),  I  =  1,  24. 

The  number  of  elements  in  this  material. 

Switch  for  type  of  equation-of-state. 

Must  be  -3  for  this  option. 

pQ  -  Reference  density. 

TQ  -  Reference  temperature. 

If  T  <  0,  code  sets  TQ  «  0.02567785  eV  (298°K). 


Variable 

(41-50) 

5. 

K,  -  Constant  in  Rosseland  opacity  expression. 

2 

Units  are  cm  /gm. 

Variable 

(51-60) 

6. 

1.  -  Integer  which  must  be  entered  but  not  used. 

Variable 

(61-70) 

7. 

a)  -  LLL  constant.  (See  variable  18  below). 

Variable 

(71-80) 

8. 

Kg  -  Constant  in  Rosseland  opacity  expression. 

Units  are  eV. 

If  Kg  <  0,  code  sets  Kg  =  0.25  eV. 

Vari able 
(1-10) 

9. 

Blank. 

Variable 

(11-20) 

10. 

Blank. 

Variable 

(21-30) 

11. 

Cy  -  Heat  capacity. 

If  Cy  <  0,  code  sets  Cy  =  3  No  k. 

Variable 

(31-40) 

12. 

Tm  -  Melting  temperature  (positive  number) 

or 

(-Em)  -  Energy  to  the  melting  point  at  zero  pressure 
from  the  reference  point  (negative  number). 

Variable 

(41-50) 

13. 

Blank. 

Variable 

(51-60) 

14. 

Blank. 

Variable  15.  H  -  Thermal  conductivity  coefficient.  If  zero, 

(61-701  0 

'  conduction  is  not  included.  Note  that  the  units 

r 

of  H  =  H  T  41  are  ergs/(cm  sec  eV). 
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I 

Variable  16. 

(71-80) 

C^i  -  Temperature  dependence  of  thermal  conduction 
coefficient  (see  variable  15). 

I 

Variable  17. 

(1-10) 

pmin 
0.8  P( 

-  Lowest  allowed  solid  density,  usually  about 
.  If  zero  or  negative,  code  sets  pm1-n  =  0.8 

i 

Variable  18. 

(11-20) 

A.  \ 

I 

Variable  19. 

(21-30) 

B.  1 

LLL  Constants. 

I 

Variable  20. 

(31-40) 

Ri  | 

(see  reference  6) 

I 

Variable  21. 
(41-50) 

R2 

I 

Variable  22. 
(51-60) 

Blank 

• 

L 

Variable  23. 
(61-70) 

Blank 

« 

1 

Variable  24. 
(71-80) 

Blank 

• 

Card  5. 


Format  (5(F5.0,  E10.3) 


There  is  one  set  of  the  following  variables  for  each  element  in  variable  1, 
card  2. 

1=1,  number  of  elements.  See  reference  7  for  element  table. 


Variable  Odd. 
Variable  Even. 


Z(I)  -  Atomic  number  of  element. 

Unnormalized  atomic  number  fraction  of  element  [C0T(I)j. 
or 

-(Unnormalized  atomic  weight  fraction  of  element.) 

All  elements  should  be  defined  in  the  same  way. 


K3  is  the  main  computation  program.  It  requires  a  variable  amount  of 
central  memory  storage  (CM)  and  extended  core  storage.  The  CM  varies  even 
for  fixed  mesh  size  since  there  are  five  sets  of  optional  subroutine  packages 
which  need  to  be  loaded  only  when  specialized  calculations  are  used.  These 
are: 

1.  elastic-plastic  package 

2.  energy-fl  ow  package 

3.  energy  source/high  explosive  package  (not  operational) 

4.  mesh  fix  package  (not  operational) 

5.  material  history  package  (not  operational) 

It  should  be  obvious  which  packages  are  required  for  a  given  computation. 

The  following  suggestions  are  included  to  aid  the  user.  Their  meaning 
will  be  clearer  after  going  through  the  input  instructions. 

1.  Each  run  by  K3  is  considered  as  a  restart  of  a  previous  run. 

K3GEN  must  be  employed  to  generate  the  first  restart  dump  for 
input  unit  1  (dump  number  1). 

2.  Cell  centers  are  tested  to  determine  if  a  cell  is  inside  a  region 
defined  as  an  energy  source. 

3.  If  the  program  immediately  has  an  address  out  of  range  error 
(MODE  =1),  check  optional  subroutine  packages  above. 

Tape  Units  K3 

1.  Restart  input  tape  and  standard  restart  output  tape  (must  be  present). 

2.  Optional  restart  output  tape. 

3.  Plot  output  tape. 

55.  Material  history  output  tape. 

67.  Internal  storage. 

77.  Internal  storage. 

Note:  Both  restart  and  plot  tapes  have  identical  formats.  The  program  can 
be  restarted  from  a  plot  tape  or  plots  can  be  made  from  a  restart  tape. 
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K3  Inputs 


Card  1 ,  Format  (7A10,  IIP) 

(1-70)  BO)  title  (MUST  agree  with  1-40  in  title  defined  in 

K3GEN,  41-70  can  be  changed). 

(71-80)  NSTART  -  Restart  tape  dump  number  to  begin  run. 

(No.  from  K3GEN  is  one)  Following  restart 
dumps  are  added  to  the  tape  after  the  last 
one  read  to  start  run.  This  can  be  changed 
by  variable  10  on  card  2.  If  NSTART  <  0 
code  will  read  to  the  last  dump  on  the 
tape  and  start  from  there.  If  NSTART  »  0 
code  sets  NSTART  *  1 . 


Card  2. 

Format  (E10.3.  1415) 

Variable 

(1-10) 

1. 

TSTOP 

*  problem  stop  time* 

If  TSTOP  =  0,  code  sets  TSTOP  =  -. 

Variable 

(11-15) 

2. 

ICYCLES 

=  problem  stop  cycle* 

If  ICYCLES  =  0,  code  sets  ICYCLES  = 

Variable 

(16-20) 

3. 

MED ITS 

=  ±  number  of  short  print  intervals  in  card  set  4 
0  <  |MEDITS|  <  9 

Variable 

(21-25) 

4. 

MED  IT 

=  ±  number  of  long  print  Intervals  in  card  set  5. 
0  <  |MEDIT|  <  9 

Variable 

(26-30) 

5. 

MOVIE 

*  ±  number  of  plot  dump  intervals  in  card  set  6. 

0  <  | MOVIE {  <  9 

Variable 

(31-35) 

6. 

NDUMPF 

*  seconds  of  C.P.  time  between  restart  tape  dumps 
If  NDUMPF  <  0  code  sets  NDUMPF  «  3600  sec. 

Variable 

(36-40) 

7. 

NDTMAX 

■  number  of  Atm,„  intervals  in  card  set  7. 
max 

0  <  NDTMAX  <  9 

*  Program  will  terminate  on  the  computer  C.P.  clock  If  TSTOP  and  ICYCLES  are 
not  defined. 


Variable  8. 
(41-45) 


Variable 

(46-50) 


Variable 

(51-55) 


Variable 

(56-60) 


Variable 

(61-65) 


Variable 

(66-70) 


Variable 

(71-75) 


Variable 

(76-80) 


Card  3. 

Variable 

(1-5) 


9. 


10. 


11. 


12. 


13. 


14. 


NDTMIN  =  number  of  At  .  intervals  in  card  set  8. 

min 

0  <  NDTMIN  <  9 

NTBOUND  =  nunber  of  radiation  boundary  temperatures 
to  be  defined  (not  operational). 

NREST  =  restart  output  switch 

If  NREST  =  0,  output  tape  is  1 

(same  as  restart  input) 

If  NREST  t  0,  output  tape  is  2. 

KNARRX  =  switch  for  card  set  9. 

If  <  0,  no  card  set  9 

If  >  0,  allows  change  of  any  variable  in 

integer  common  block  with  card  set  9. 

KARRX  =  switch  for  card  set  10. 

same  as  KNARRX  except  for  floating  common 
block  and  card  set  10. 

NEXT VI  =  print  switch  for  nonstandard  edit  variable.. 

See  card  4  in  K3GEN  inputs  but  here  one  inputs 
the  value  of  the  zone  storage  order  variable, 
i.e.,  value  of  LK(NEXTVl)  as  in  K3GEN.  Use 
only  if  it  is  desired  to  change  K3GEN  inputs. 

!€XTV2  *  same  as  NEXTV1  for  second  printed  variable. 


15. 


JEXTV3  =  same  as  NEXTV1  for  third  printed  variable. 


Format  (4F5.0.  E10.3,  9F5.0) 

1.  SOURCES  =  switch  for  energy  sources 

5  0  no  sources 

*  1  input  source  data  (not  operational) 
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Variable  2. 
(5-10) 


Variable  3. 
(11-15) 

Variable  4. 
(16-20) 


Variable  5. 
(21-30) 


Variable  6. 
(31-35) 


Variable  7. 
(36-40) 


Variable  8. 
(41-45) 


ELPL  *  switch  for  elastic-plastic  calculation. 

■  0  to  turn  off  EP  calculation. 

>  1  EP  calculation  on  and  card  set  12  is  required. 
The  value  of  ELPL  determines  how  yield  functions 
are  treated  in  a  mixed  cell. 

*  1,  volume  average  yield  parameters  (normal) 

=  2,  =  0  in  mixed  cell. 

*  3,  Y  *  0  in  mixed  cell . 

POROUS  Not  Used 

FRACT  =  switch  for  input  of  fracture  calculation  data. 

See  card  set  11.  If  FRACT  ?  0  input  card  set  11. 
If  FRACT  =  1  input  pressures. 

If  FRACT  =  2  input  stresses. 

GRAV  ■  acceleration  of  gravity  (980  cm/sec)  if  desired; 

zero  otherwise,  (acceleration  in  -  Y  direction 
only). 

PR  INTO  =  switch  to  control  long  edit  on  1st  cycle  of 

current  run 

=  0  long  edit  on  1st  cycle 
j*  0  no  long  print  of  1st  cycle  (but  gives  short 
edit). 

CPSHIFT  *  variable  to  allow  for  extra  CP  time  in  seconds 
in  computer  clock  run  termination.  Used  to 
allow  time  for  tape  copying,  plots,  etc.  Default 
value  Is  15. 

FNTRACE  *  switch  to  allow  changes  In  tracer  particle 
position 
>0.  no  change 

i  0.  input  data  (not  operational). 


--W/ 


mipm  -  ■***  •»»* 
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Variable  9. 
(46-50) 


Variable  10. 
(51-55) 


Variable  11 . 
(56-60) 


Variable  12. 
(61-65) 


Variable  13. 
(66-70) 


Variable  14. 
(71-75) 


EXTAPM  =  switch  to  allow  continuation  of  plot  output 
file  on  previously  generated  reel  (unit  3). 

See  variable  5,  card  2  and  card  6. 

=  0,  generate  new  file  on  output  unit  3 

>  0,  add  new  data  after  data  generated  during 
a  previous  run.  EXTAPM  is  the  last  cycle 
number  to  be  retained  from  old  data.  Cycles 
with  cycle  number  >  EXTAPM  are  overwritten. 

<  0.  add  new  data  after  last  cycle  dump  on  the 
old  file.  (EXTAPM  *  0  for  cycle  0  always) 

FRACDT  -  Fraction  of  Courant  stability  used  to  calcu- 

Atcs*  Nomally  °*6»  no  case  9reater  than 
one.  If  FRACDT  _<  0,  code  sets  FRACDT  *  0.6. 

DTI NCR  -  Maximum  factor  allowed  for  increase  in  at 
from  one  cycle  to  the  next. 

If  DTINCR  <  0,  code  sets  DTI NCR  *  1.068000433 

FIXMH  -  Switch  to  allow  definition  of  fixed  mesh  lines. 

=  0,  not  used. 

j4  0,  input  data  (not  operational) 

HIST  *  Switch  for  output  on  history  data  tape  (unit  55). 

»  0,  no  output 

f  0,  output  data  to  tape  55.  Output  frequency 
is  every  )HIST|—  cycle.  (Not  operational) 

RH0ED  -  Switch  for  extended  long  edit  giving  density 
and  volume  fraction  for  each  material. 

*  0.  no  extended  edit 

*  1.  produce  extended  edit  (used  only  if  more 
than  two  materials  present.) 
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Card  Set  4  Short  Print  Interval 

Use  either  4.1  or  4.2  depending  on  MEDITS  on  card  2.  There  are  |MEDITS| 

sets  of  data  pairs. 

If  MEDITS  >  0 

Card  4.1  Format  (8E10.3) 

I  -  1 ,  MEDITS 

Variable  Odd  -  PTIMES(I)  3  time  at  which  the  short  edit  frequency 

changes  from  DPTIMES(I-l)  to  DPTIMES(I). 
[PTIMES(l)  .<  start  time  of  problem] 

Variable  Even  -  DPTIMES(I)  3  the  time  interval  between  short  edits 

from  PTIMES(I)  to  PTIMES(I+1) 

For  times  >  PTIMES(MEDITS) ,  interval  UPTIMES (MED ITS)  is  used  to  the  end  of 
the  calculation. 

If  MEDITS  <  0 

Card  4.2  Format  (8110) 

Same  as  4.1  except  cycles  are  used  instead  of  time.  The  variables  are  named 
NCYEDS(I)  and  NDCYEDS(I). 


Card  Set  5  Long  Print  Interval 

Use  either  5.1  or  5.2  depending  on  MEDIT  on  card  2.  There  are  |MEDIT|  sets 

of  data  pairs.  A  short  edit  is  always  given  as  part  of  the  long  edit. 


Warning — long  edits  are  long  for  large  problems. 


If  MEDIT  >  0 

Card  5.1  Format  (8E10.3) 

Same  as  card  4.1  except  for  long  edit.  The  variable  names  are  PTIME(I) 
and  DPTIME(I). 


If  MEDIT  <  0 

Card  5.2  Format  (8110) 

Same  as  card  4.2  except  for  long  edit.  The  variable  names  are  NCYED(I) 
and  NDCYED(I). 

Card  Set  6  Plot  Tape  Dump  Intervals 

Present  only  if  MOVIE  t  0  on  card  2. 

Use  either  6.1  or  6.2  depending  on  MOVIE  on  card  2.  There  are  (MOVIE |  sets 
of  data  pairs. 

If  MOVIE  >  0 

Card  6.1  Format  (8E10.3) 

Same  as  card  4.1  except  for  plot  tape  dumps.  The  variable  names  are  TMOVIE(I) 
and  DTMOVIE ( I ) . 

If  MOVIE  <  0 

Card  6.2  Format  (8110) 

Same  as  card  4.2  except  for  plot  tape  dumps.  The  variable  names  are  IM0V1E(I) 
and  IDMOVIE (I ) . 

Card  Set  7  Format  (8E10.3)  At  .  data 

max 

Present  only  If  NDTMAX  >  0  on  card  2. 

Same  form  as  card  4.1.  The  variable  names  are  TDTMAXX(I)  and  DTMAXX(I). 

Card  Set  8  Format  (8E10.3)  Atm1n  data 

Present  only  If  NDTMIN  >  0  on  card  2. 

Same  form  as  card  4.2.  The  variable  names  are  TDTMINN(I)  and  DTMINN(I). 


Card  Set  9  Format  (2110) 

Present  only  If  KNARRX  >  0  on  card  2. 

This  card  set  allows  the  changing  of  a  variable(s)  in  the  integer  COMMON 
block.*  CAREFUL 

Variable  1.  I  *  location  in  NARRX  array. 

(1-10) 

Variable  2.  NARRX(I)  »  new  value  to  put  in  location  NARRX(I) 

(11-20) 

Repeat  cards  until  finished. 

Insert  blank  card  to  stop  reading  these  cards. 


Card  Set  10.  Format  (110,  E10.3) 

Present  only  If  KARRX  >  0  on  card  2. 

This  card  set  allows  the  changing  of  a  variable(s)  in  the  floating  COMMON 
block.*  CAREFUL 

Variable  1.  I  *  location  In  ARRX  array. 

(1-10) 

Variable  2.  ARRX(I)  »  new  value  of  ARRX(I). 

(11-20) 

Repeat  cards  until  finished.  Insert  blank  card  to  stop  reading  these  cards. 


♦For  example  IACT  -  NARRX(17) 
TIME  -  ARRX (33) 
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Card  Set  1 1 . 


Fracture  Data 


Present  only  If  FRACT  +  0  on  card  3. 

If  FRACT  -  1,  all  Inputs  are  pressure  values. 

If  FRACT  «  2,  all  Inputs  are  stress  values  except  for  P, 
pressure. 

values  []  are  used  for  any  zero  input  value  or  if  FRACT  =  0.  Since  positive 
numbers  have  no  meaning  in  the  present  model ,  the  code  will  add  the  minus  if 
a  positive  number  is  found. 


fRACT  VOID  wh1ch  1S 

All  quantities  are  negative  and  in  units  of  dynes/cm2.  Default 


Card  11.1 
Variables 
K  *  1 . NUMMT 


Variable 


Variable 


Variable 


Format  (8E10.3) 


(K) 


TRACT 


TRACT  MIX 


TRACT  VOID 
ACES 


minimum  stress  or  pressure  in  a  cell  of 

material  K  with  no  void  present.  There 

is  one  input  for  each  material.  [-2xl010] 

minimum  stress  or  pressure  in  a  cell  of 

mixed  material  with  no  void  present. 

Inputs  on  card  11.2  will  override  this 

8 

variable  for  special  mixtures.  [-1x10  ] 
minimum  pressure  in  any  cell  containing 
void.  [-5x10®] 


*  number  of  cards  11.2  to  define  interface 
strengths.  An  integer  but  punched  with 
decimal  point  0.  <  NFACES  ±  12.  (floating 
number) 

Continue  on  card  with  identical  format  for  6  or  more  materials. 


Card  11.2  Format  (215,  E10.3) 

There  are  NFACES  (card  11.1)  of  these  cards  to  define  strength  at  material 
interfaces  between  any  two  materials.  I  *  1,  NFACES 


Variable 

1. 

MATF1 

«  one  of  materials  at  Interface 

(1-5) 

Variable 

2. 

MATF2 

■  other  material  at  interface 

(6-10) 

Variable 

(11-20) 

3. 

PFRACT 

■  minimum  stress  or  pressure  In  a  cell 
containing  only  materials  MATF1  and 
MATF2.  (no  void  or  other  materials) 

tPfRACT  MIX3 
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Card  Set  12.  Format  (7E10.3)  Elastic  Properties 


Present  only  If  ELPL  >  0  on  card  3. 

The  elastic- plastic  package  of  sub  routines  must  be  loaded  with  K3  if  this 
option  is  used. 

SHOULD  be  included  for  EP  calculation 

one  of  the  following  cards  for  each  of  the  NUMMT  materials.  Material 
etc.  Input  blank  card  for  hydrodynamic  material. 

1.  Yq  (dynes/cm^) 

2.  em  =  absolute  value  of  melt  energy  (ergs/gm) 

If  not  defined,  the  program  will  determine  for  EOS. 


There  is 
1  first. 

Variable 

(1-10) 

Variable 

(11-20) 


Variable  3. 
(21-30) 

Variable  4. 
(31-40) 


vQ  =  Poisson  ratio. 

Csmin  =  nrin''mum  sound  velocity  to  be  used  to 
calculate  shear  modulus,  (normally  zero) 


Variable  5. 
(41-50) 


Variable  6. 
(51-60) 


a  *  fraction  of  melt  energy  where  material  begins  to 
lose  strength. 

If  a  <  0,  code  sets  a  *  0.8 

RHOEPU  »  py.  If  Pj  <.  0,  code  sets  pu  =  0.9  pm,  where  pm 
is  the  liquid  density  at  zero  pressure  melt. 


Variable  7. 


RHOEPL  =  pL<  If  pL  <  0,  coee  sets  p^  *  0.9  pL. 


Card  Set  13.  Reserved  for  future  expansion. 


Uppw  |,W«  I  y 
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4.  K3PLT 


a.  General 

K3PLT  is  a  program  used  to  produce  dot  plots  on  the  U.S.  Army's  Advanced 
Research  Center's  interactive  color  graphics  system.  (See  Reference  5  for  a 
description  of  the  system  and  operating  instructions).  K3PLT  is  a  modifica¬ 
tion  of  CSQPLT  and  the  SCORS  support  plotting  package  and  produces  a  two- 
dimensional  image  of  selected  parameters  edited  from  K3  restart/ dump  tapes. 

The  main  use  of  K3PLT  is  to  produce  material  density  plots  but  can  also, 
without  modification,  plot  pressure  or  temperature.  A  subroutine  OPUS  is 
provided  allowing  the  user  to  modify  the  program  with  minimal  effort  in  order 
to  select  other  parameters  for  plotting. 

The  amount  of  extended  core  storage  (ECS)  currently  programmed  into  K3PLT 
for  storage  of  plot  variables  in  100010.  This  amount  can  be  increased  or  re¬ 
duced,  if  desired,  to  suit  problem  requirements  and  thereby  increasing  program 
computer  efficiency.  The  amount  of  ECS  required  is  4*IMAX*JMAX*KMAX,  but  can 
get  by  with  less,  using  disk  file  for  excess  storage  requirements. 

As  stated  above,  K3PLT  is  an  outgrowth  of  the  CSQPLT  program.  The 
modifications  made  have  negated  some  of  the  capabilities  of  the  CSQPLT  pro¬ 
gram  and  as  a  result,  some  Inputs  necessary  for  CSQPLT  are  no  longer  required 
for  K3PLT.  These  inputs  are  included  in  the  description  of  K3PLT  for  com¬ 
pleteness  sake,  and  in  the  event  at  some  future  time  it  is  deemed  useful  to 
Include  the  former  capabilities. 

b.  Input  Instructions 

K3PLT  Is  run  Interactively,  with  the  user  inputting  control  parameters 
via  the  system  keyboard  and  trackball  control.  The  following  instructions 
lead  the  user  step-by-step  through  the  procedures  necessary  for  obtaining 
pictures  (plots)  on  the  system's  CRT. 

Step  1  After  job  submittal  a  logo  will  appear  on  the  screen  of  an  available 
station,  with  a  blinking  cursor  (the  logo  Is  currently  "KAMAN  SCIENCES  CORP¬ 
ORATION"  on  a  blue  background).  The  cursor  should  be  positioned  onto  the 
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lower  quarter  of  the  screen*  using  the  trackball,  and  the  trackball  ENTER 
button  depressed.  This  will  advance  the  program  to  the  next  step,  definition 
of  inputs. 

Step  2  The  screen  will  display  "ENTER  OATA  FOR  CARD  NO.  1"  and  "K0N(**)=". 
The  user  should  then  enter  the  variable  number  (see  Figure  IV-1,  list  of 
card  1  inputs),  via  the  keyboard,  of  KON  that  is  to  be  changed  from  the 
default  value.  These  inputs  can  be  in  any  order  desired.  After  the  variable 
number  is  selected  (which  must  be  an  integer),  the  screen  will  then  display 
"K0N(n* )=*♦********".  The  value  of  KON(n)  should  then  be  entered  (which  also 
must  be  an  integer).  The  program  will  cease  expecting  values  of  KON  when 
K0N(16)  is  defined,  or  the  user  may  skip  to  the  next  step  by  entering  a  slash 
(/)  for  any  value  of  KON*^.  If  the  program  is  to  be  terminated,  "ENDRUN" 
should  be  entered  for  any  KON.  Once  the  inputs  are  completed,  the  cursor 
will  then  blink,  and  the  user  may  advance  to  the  next  step  by  depressing  the 
trackball  ENTER  button*. 


t  When  the  cursor  is  positioned  in  the  upper  three-quarters  of  the  screen 
and  the  ENTER  button  depressed,  the  cursor  position  identifies  the  location 
at  which  alphanumeric  data  Is  to  be  entered  from  the  keyboard  (1  to  10  charac¬ 
ters).  This  Is  handy  for  annotating  plots  for  future  identification.  After 
entering  the  data  (depressing  the  keyboard  ENTER  button),  position  the  cursor 
to  the  lower  portion  of  the  screen  and  proceed. 

ft  For  values  of  KON  not  defined,  either  the  default  values  will  be  used, 
or  values  of  KON  previously  defined. 
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KON 

Default  Value 

Used 

Description 

1 

1 

Yes 

*  Input  tape  number  (1-5) 

2 

0 

No 

=  0,  display  on  CRT 
=  1 ,  user  tape  40  output 

3 

3 

Yes 

=  0,  velocity  vector  (not  available) 

*  1,  momentum  vector  (not  available) 

=  2,  pressure 
=  3,  density 

«  4,  user  option  (see  subroutine  OPUS- 
temperature  will  be  plotted  as 
presently  coded) 

4 

1 

Yes 

=  0,  plot  all  tape  dumps 
i  0,  plot  only  selected  dumps  (see  also 
variables  13  and  14) 

5 

1 

Yes 

=  number  of  mesh  frames  and  number  of  dot 
scale  frames  (>.0)  -  automatically  reduced 
to  zero  after  first  plot 

6 

1 

Yes 

=  number  of  title  frames  (_>  0) 

7 

1 

Yes 

=  0,  no  mirror  image  (in  X  direction) 

=  1,  produce  mirror  image 

8 

1 

Yes 

=  number  of  times  to  produce  first  plot  (>^  1) 

9 

1 

Yes 

=  number  of  times  to  produce  following  plots 
(>  1) 

10 

0 

No 

=  switch  to  plot  interface  1  ines/ symbols 
(not  available  -  should  be  zero) 

Figure  IV-1 

K3PLT  Inputs  -  Card  Set  1 
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ICON 

Default  Value 

Used 

Description 

11 

1 

Yes 

=  number  of  THE  END  frames  (>.  0) 

12 

0 

No 

*  switch  for  frame  size  (not  available  - 

should  be  zero) 

13* 

1 

Yes 

=  first  tape  dump  to  plot  (see  variable  4) 

14* 

20000 

Yes 

-  last  tape  dump  to  plot  (see  variable  4) 

15 

0 

No 

=  0,  produce  dot  plots 

>.  1,  produce  contour  plots  (not  available) 

16 

0 

No 

3  0,  do  not  plot  tracer  particles 

=  1,  plot  tracers  with  other  plot  data 
(not  available) 

=  -1,  plot  tracers  without  other  plot  data 
(not  available) 


*  Program  uses  these  inputs  even  if  K0N(4)  is  zero. 


r 


\ 
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Figure  IV-1  (continued) 
K3PLT  Inputs  -  Card  Set  1 
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Step  3  The  program  will  ask  if  the  user  wishes  to  define  two  titles  for 
display  on  the  title  frame  to  be  produced  at  a  later  stage  of  program  exe¬ 
cution.  “YES"  should  be  entered  If  the  user  wishes  to  enter  titles  on  the 
first  title  frame.  Any  other  entry  will  produce  an  unlabeled  title  frame 
(first  title  frame  only)  and  advance  the  program  to  the  next  step.  If  titles 
are  to  be  entered,  the  keyboard  is  used  to  enter  the  desired  alphanumeric  data. 
Again,  once  the  program  digests  the  user  inputs,  a  blinking  cursor  signals  that 
the  program  Is  waiting  to  be  advanced. 

Step  4  The  screen  will  then  display  "ENTER  DATA  FOR  CARD  NO.  3"  followed  by 
XBL  a  **********".  Refer  to  Figure  IV-2  for  the  definition  of  these  input 
variables.  Entering  a  slash  will  automatically  set  the  default  values,  or  the 
values  last  used,  for  the  displayed  variable  and  any  variables  remaining  to  be 
defined  for  this  input  set.  The  values  entered  must  be  in  floating  print  for¬ 
mat,  with  or  without  an  Integer  exponent.  The  following  are  examples  of  allowed 
Inputs: 

1.0 

2.E-2 

1.99999 

-6.89E20 


After  all  variables  are  defined,  or  a  slash  entered,  the  cursor  will  blink, 
and  the  trackball  ENTER  button  should  be  depressed  to  advance  to  the  next 
step. 


Variable 

Default 

Description 

XBL 

0.0 

X  -  left,  (left-most  value  of  abscissa)  - 
automatically  set  to  minus  X  -  right  if 

K0N(7)  is  equal  to  1. 

XBR 

0.0 

X  -  right 

YBB 

0.0 

Y  -  bottom 

YBA 

0.0 

Y  -  top 

ZBB 

0.0 

Z  -  back.  Values  of  K-planes  less  than 
this  will  not  be  displayed 

ZBF 

0.0 

Z  -  front.  Values  of  K-planes  greater 
than  this  will  not  be  displayed 

ROTA(l)* 

0.0 

Angle  (in  degrees)  of  X,  Y,  and  Z  rotation 

R0TA(2)* 

0.0 

of  density  plots.  X  rotation  is  performed 

R0TA(3)* 

0.0 

first  followed  by  Y,  then  Z.  Rotation  is 
about  the  midpoint  of  X  -  left  and  X  -  right, 

Y  -  bottom  and  Y  -  top,  Z  -  back  and  Z  -  front 

SCALEM 

0.0  | 

For  velocity  and  momentum  plots  only  -  not 

VECTMIN 

0.0  | 

avail  able 

DOTFM 

2.0  ] 

Dot  density  control,  given  by 

DOTPOW 

1.0  j 

[  Variable]  D0TP0W 
|_  “DOTFM  J 

DOTPOW  is  restricted  to  0.<  DOTPOW  <  1.0 

*  if  rotations 

are  used,  grid  labeling  is  meaningless 

Figure  IV-2 

K3PLT  Inputs  -  Card  Set  3 


Step  5  If  selected  dumps  are  to  be  plotted  (see  definition  of  K0N(4)),  the 
program  will  expect  the  user  to  define  the  first  dump  to  be  plotted.  Also, 
after  the  completion  of  each  plot,  the  program  will  return  to  this  step  (if 
K0N(4)  is  not  zero)  at  which  time  the  next  dump  to  be  plotted  may  be  selected. 
If  a  zero  is  entered,  the  program  will  return  to  step  2,  at  which  time  the 
run  may  be  terminated  or  new  values  of  selected  inputs  may  be  entered.  Note 
that  dumps  selected  must  be  in  ascending  order  -  if  a  dump  of  lower  numerical 
value  is  to  be  plotted,  a  return  to  step  2  must  be  made. 

After  the  dump  number  is  entered,  the  cursor  will  blink,  and  after  the 
ENTER  button  is  depressed,  the  program  will  initiate  the  actual  plotting 
sequence. 

c.  PI  ots 

The  program  will  display  a  title  frame,  a  mesh  frame,  and  a  plot  of 
dot  densities  (if  the  appropriate  KON  values  are  non-zero)  before  the  first 
dump  is  plotted.  After  each  frame  is  displayed  the  cursor  will  blink  sig¬ 
nalling  the  frame  is  complete  and  the  program  is  waiting  before  advancing 
to  the  next  frame.  This  allows  the  user  time  to  make  hardcopy  prints  of 
the  display  or  to  study  the  screen.  Advancement  to  the  next  frame  is  made 
by  depressing  the  trackball  ENTER  button  (see  also  the  footnote  with  step  1). 

The  displays  of  selected  dumps  are  built  by  displaying  one  K-plane  at 
a  time.  The  trackball  cursor  will  blink  after  each  plane  is  displayed.  The 
ENTER  button  must  be  depressed  to  display  the  next  K-plane.  If  the  user 
wishes  to  terminate  the  display  (not  display  any  more  K-planes  for  this  plot) 
the  trackball  cursor  can  be  positioned  in  the  upper  three-quarters  of  the 
screen,  the  ENTER  button  depressed,  then  the  trackball  cursor  returned  to 
the  lower  portion  of  the  screen  (where  it  will  again  blink,  the  program 
waiting  to  advance).  After  all  K-planes  have  been  displayed,  or  the  display 
terminated  by  proper  positioning  of  the  trackball  cursor,  the  program  will 
return  to  step  5  If  K0N(4)  is  not  zero  or  display  the  next  dump.  After  all 
dumps  have  been  displayed,  a  "THE  END"  frame  will  be  displayed  (if  KON(ll) 
is  non-zero),  and  the  program  will  return  to  step  2. 
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APPENDIX  A 
K3  TEST  RUNS 


The  purpose  of  this  section  is  to  provide  the  K3  user  with 

o  an  overview  of  the  runs  made  to  test  and  debug  the 
K3  family  of  programs, 

o  example  inputs  for  a  specific  problem, 

o  example  outputs  from  the  K3PLT  program 

As  can  be  imagined  by  those  familiar  with  large  program  development, 
many  unsuccessful  test  runs  were  made  before  satisfactory  results  were 
obtained.  Plus,  the  test  problems  were  set  up,  not  necessarily  to  model 
a  physically  realistic  situation,  but  rather  to  test  out  a  specific  por¬ 
tion  of  the  code.  These  runs  are  not  described  here.  What  is  described 
include,  first,  a  series  of  runs  made  to  test  a)  the  symmetry  of  results, 
i.e.,  does  the  code  predict  material  response  and  flow  in  a  given  direction 
consistent  with  that  for  the  other  directions,  and  b)  that  results  are 
similar  to  those  predicted  by  CSQ  for  the  two-dimensional  analog.  These 
test  runs  were  configured  with  a  27  cm3  cube  (three  centimeters  to  a  side) 
of  high  density  gas  situated  in  the  center  of  a  7x7x7  centimeter  box.  The 
motion  of  the  gas  was  observed  as  the  gas  expanded  to  fill  the  container. 

Secondly,  there  is  described  a  series  of  runs  that  were  made  with  a 
physically  realistic  problem  (indeed,  the  problem  models  an  actual  physical 
experiment)  in  order  to  compare  results  with  those  predicted  with  CSQ,  to 
test  out  a  large  problem,  and  to  test  the  completely  assembled  code.  Since 
this  last  problem  is  of  more  Interest  than  others,  it  is  used  to  cover  the 
last  two  items  above,  i.e.,  example  Inputs  and  outputs.  The  problem  was  to 
simulate  the  response  of  a  block  of  aluminum  with  a  steel  backing  plate 
Impacted  by  a  steel  ball  moving  approximately  5000  meters  per  second. 

1 .  Gas  in  a  Box  Problem 

The  Initial  configuration  of  the  problem  is  shown  in  Figure  A-l, 
where  a  3x3x3  grid  of  gaseous  material  is  placed  in  the  center  of  a 
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7x7x7  grid  box.  Initially  the  sides  of  the  box  were  solid  to  test  boundary 
conditions.  The  reasons  for  choosing  the  above  grid  sizes  were: 

o  To  reduce  run  time  and  amount  of  output  by  using  the  minimum 
size  grid, 

o  The  minimum  size  of  any  active  region  should  be  at  least 
3  cells  in  any  direction, 

o  To  provide  space  for  the  gas  to  expand,  of  a  size  to  reduce 
computation  time  yet  allow  sufficient  space  to  study  movement 
from  cell-to-cell  before  boundary  contact, 

o  Since  CSQ  in  the  Cartesian  mode  assumes  a  one  centimeter  thick 
cell  in  the  imaginary  z  direction,  one  centimeter  spacing  for 
the  entire  grid  was  chosen  to  retain  uniformity,  for  ease  in 
interpreting  output,  and  for  direct  comparison  purposes. 

The  gas  was  polyethelene  (this  choice  was  arbitrary  -  any  gas  would  have 

been  sufficient  for  the  problem)  with  a  density  of  2  grams/cm  and  an 
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internal  pressure  of  5  x  10  dynes/cm.  The  total  mass  was  54  grams. 

With  the  void  fraction  option  turned  off,  the  gas  would  reach  the 
walls  of  the  container  within  five  cycles  of  a  K3  run,  or  within  .04  micro¬ 
seconds  of  problem  time.  The  CP  time  for  a  K3  run  to  this  point  was  approxi¬ 
mately  three  seconds  which  greatly  expedited  the  turn-around  time  for  a  job. 
Results  were  deemed  satisfactory  when  the  following  was  observed: 

o  Total  momentum  of  the  system  to  the  time  of  boundary  contact 
was  zero; 

o  Velocities,  temperatures,  pressures,  etc.,  in  cells  of 

symmetrically  opposite  positions  were  of  the  same  absolute 
value  and  with  the  correct  sign; 

o  Mass  was  conserved,  i.e.,  mass  did  not  leak  through  the 
reflecting  walls  of  the  container; 


o  Total  energy  of  the  system  was  conserved  to  a  few  percen¬ 
tage  points  error; 
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o  Values  for  cell  parameters  in  the  K  =  4  plane  agreed  with 
a  CSQ  benchmark  run  throughout  the  history  of  a  run. 

Variations  on  this  configuration  were  run  to  test  the  reflective-trans- 
mittive  boundary  options,  i.e.,  some  walls  were  reflective  while  others 
were  transmitti ve,  allowing  mass  to  flow  out  of  the  problem.  Comparisons 
were  made  with  CSQ  benchmark  runs  and  with  previous  K3  runs  to  determine 
that  when  opposite  walls  were  transmittive  or  reflective,  symmetry  was 
preserved  and  that  results  were  independent  of  which  opposite  boundarys 
were  transmittive. 

Another  variation  of  this  configuration  was  formulated  to  test  out 
the  elastic-plastic  package  of  routines.  In  this  configuration,  the  gas 
was  replaced  with  a  solid  block  of  aluminum  and  given  an  initial  velocity 
in  first  the  minus-X  direction,  then,  in  another  run,  in  the  plus-Y  direc¬ 
tion,  etc.,  and  results  of  the  impact  with  the  walls  of  the  problem  examined. 
Again,  results  were  satisfactory  when  symmetry  was  preserved  and  when  they 
agreed  with  CSQ  benchmark  runs. 

2.  Aluminum  Block  -  Steel  Ball  Impact  Problem 

This  problem  was  an  attempt  to  duplicate  a  previous  CSQ  calculation 
where  CSQ  was  run  to  simulate  an  experiment  conducted  at  Arnold  Engineering 
Development  Lab,  Tulahoma,  Tennessee.  Since  the  CSQ  run  gave  excellent 
results,  in  complete  agreement  with  empirical  data,  it  was  expected  that 
K3  should  give  comparable  results  before  it  was  certified  operational. 

The  following  paragraphs  present  the  history  of  the  K3  run  and  comparisons 
with  the  original  CSQ  run. 

The  problem  as  initially  configured  is  shown  in  Figure  A -2a  for  CSQ 
and  Figure  A-2b  for  K3  (neither  figure  is  to  scale).  Units  are  in  centi¬ 
meters.  CSQ  was  run  using  the  cylindrical  geometry  option,  therefore  the 
aluminum  block  (and  steel  backing)  as  well  as  the  steel  ball  were  taken 
as  figures  of  revolution  about  the  Y-axis  with  radii  of  15.24  and  .9166 
centimeters,  respectively.  K3  was  run  using  the  only  available  option, 
Cartesian  geometry,  with  reflective  boundarys  on  the  X  =  0  and  Z  =  0 
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planes  (Y  -  Z  plane  and  X  -  Y  plane).  This  essentially  reduced  the  problem 
configuration  to  one- fourth  of  the  physical  geometry,  i.e.,  only  that  part 
of  the  block  and  sphere  which  is  in  the  positive-X,  positive-Z  quadrant  is 
simulated.  The  dimensions  for  the  K3  configuration  were  identical  to  those 
for  CSQ. 

For  the  CSQ  run,  there  were  NVAR  =  16  variables  stored  for  each  cell 
in  the  problem*,  or  a  total  of  I*J*NVAR  =  92960  words.  K3  was  zoned 

ilia  a  lildX 

with  ^ax  =  18,  Jmax  =  35»  and  Kmax  =  18»  and  with  NVAR  =  17»  used  a  tota1 
of  192,780  words  of  storage  for  cell  variables.  This  particular  zoning 

kept  the  required  size  of  cell  storage  below  the  size  which,  if  exceeded, 

would  force  some  storage  onto  disk  and  increase  the  run  time.  Though  the 

zoning  was  more  coarse  for  the  K3  problem,  previous  studies  using  CSQ  with 

fine  and  coarse  grids  showed  the  size  of  the  zones  had  little  effect  on  the 

results  for  problems  of  this  type,  as  long  as  the  spacing  was  judiciously 

chosen. 

The  K3GEN  inputs  for  this  problem  are  shown  in  Figure  A-3a.  Note  the 
velocity  of  the  ball  is  4.995  x  105  cm/sec  in  the  +Y  direction.  The  K3 
inputs  are  shown  in  Figure  A-3b.  K3  was  set  to  run  200  cycles,  which  trans¬ 
lated  to  approximately  16  ys  (computer  time  required  was  slightly  over  one- 
half  hour  of  CP  time).  The  restart  tape  produced  by  K3  was  then  input  to 
K3PLT  and  selected  dumps  plotted.  Figures  A-4,  A-5  and  A-6  show  the  title, 
mesh  and  dot  scale  frames  produced  (by  user  options)  by  K3PLT.  As  the  mesh 
frame  indicates,  the  reflective  option  in  K3PLT  was  chosen,  and  the  plot 
scale  chosen  so  -8.0  j<  x  8.0  i.e.,  only  a  portion  of  the  aluminum  block/ 
steel  backing  is  shown,  effectively  enlarging  the  area  of  interest.  Figure 
A-7  shows  the  configuration  of  the  K3  problem  at  time  zero,  K  =  1  plane  only. 


*  Three  materials  were  defined,  since  the  steel  backing  and  steel  ball 
were  different  alloys.  In  K3,  in  order  to  reduce  NVAR,  the  steels  were 
assumed  to  be  alike. 
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K3GEN  INPUTS 
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Figure  A- 5 
K3PLT  Mesh  Frame 
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Figure  A- 7 

Material  Density  Plot 
at  Time  =  0.0 
K  =  1  Plane  Only 
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Figure  A-8  shows  the  K  =  1  plane  at  problem  time  4.5  us,  and  Figure  A-9 
shows  the  K  =  1  through  5  planes  at  the  same  time  (the  K-planes  overlay, 
or  are  superimposed  upon,  previously  plotted  K-planes,  hence  the  darker 
the  plots,  showing  increased  density).  Material  density  for  cycle  136, 
problem  time  7.5  ps,  is  shown  in  Figures  A-10  through  A-12.  Figure  A-10 
shows  the  K  =  1  plane,  Figure  A-ll  K-planes  1  through  8,  and  Figure  A-12 
K-planes  one  through  15.  Figures  A-13  and  A-14  show  the  material  density 
for  the  K  =  1  plane  and  K  =  1  through  8  planes,  respectively,  near  the  end 
of  the  problem.  Figures  A- 1 5  and  A- 1 6  show  a  different  view  of  K3  results 
for  this  problem  time.  Here,  only  the  right-half  of  the  configuration  is 
shown  (positive  X,  2  quadrant)  rotated  -30  degrees  about  the  X-axis  and 
+45  degrees  about  the  Y-axis.  Both  show  an  increasing  number  of  K-planes. 
Figure  A-17  gives  a  schematic  representation  of  these  views  for  better 
interpretation.  Note  the  area  of  impact  is  on  the  back-left  corner  of 
the  material  block. 

The  size  of  the  cavity  predicted  by  CSQ  and  K3  at  approximately  10  ps 
is  shown  in  Figure  A-18.  Considering  there  is  a  difference  in  time  of  1.14 
u s  between  the  CSQ  plot  and  the  K3  plot,  and  there  is  still  appreciable 
movement  in  the  mass,  the  shapes  of  the  cavities  are  in  good  agreement.* 

The  lip  formed  on  the  cavity  boundary,  however,  is  somewhat  smaller  in  K3 
than  it  is  in  CSQ.  An  examination  of  the  total  output  from  K3,  i.e.,  an 
examination  of  K-planes  away  from  the  K  =  1  plane  revealed  that  the  lip 
is  larger  as  one  progresses  away  from  the  reflective  boundary,  and  is  in 
agreement  with  CSQ  results.  This  result  may  be  due  to  conditions  imposed 
by  the  boundaries  and  the  problem  configuration  (size  of  zones),  but  fur¬ 
ther  explanation  must  await  a  detailed  study.  The  mass  in  the  extremes  of 
the  lip  is  very  small  though,  representing  fractured  or  vaporized  material 
no  longer  connected  to  the  bulk  of  the  aluminum  block,  and  therefore  has 
little  if  any  effect  on  this  particular  problem. 


*  The  plotted  points  are  shown  with  error  tolerances  which  indicate  the 
width  of  a  cell  at  the  point  and  the  uncertainty  of  locating  the  exact  mass 
boundary.  Note  the  K3  cell  widths  are  much  greater  than  the  CSQ  cell  widths 
due  to  the  coarseness  of  the  K3  grid. 
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Figure  A-10 
Material  Density  Plot 
at  Time  =  7.5  ys 
K  =  1  Plane  Only 
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Figure  A-ll 
Material  Density  Plot 
at  Time  *  7.5  ys 
K  =  1  through  8  Planes 
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Figure  A-13 
Material  Density  Plot 
at  Time  *  9.8  ys 


FIGURE  A-17 

SCHEMATIC  REPRESENTATION  OF  FIGURES  A-15  AND  A-16 
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Figure  A-19  shows  the  comparative  cavity  sizes  at  or  near  16  us 
problem  time.  With  a  close  match  in  problem  times  between  K3  and  CSQ, 
the  cavity  sizes  are  also  in  close  agreement.  Though  the  cavity  lips 
are  more  closely  matched  than  in  Figure  A-18,  there  is  a  different  shape. 
Again,  this  may  be  due  to  the  coarseness  of  the  grid. 

While  the  CSO  code  was  run  out  to  almost  100  ps,  K3  was  terminated 
at  this  point  since  little  information  could  be  gained  by  running  it  fur¬ 
ther  with  this  particular  zoning.  Future  work  on  the  code  will  be  directed 
at  incorporating  the  rezoning  feature  now  in  CSQ  into  K3,  thereby  allowing 
problems  of  this  type  to  be  run  further  with  a  reduced  execution  time. 
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